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CHAPTER  I 


INTRODUCTION  AND  SUMMARY 

The  research  undertaken  in  this  program  addresses  two  new  approaches  for 
processing  ceramic  materials.  One  is  based  on  particulate  processes  using 
powders  with  an  unusual  and  highly  specific  set  of  characteristics.  The 
other  attempts  to  modify  material  properties  by  means  of  rapid  solidification 
techniques. 

It  was  hypothesized1  and  now  has  been  demonstrated2  that  the 
micros  true tural  characteristics  of  polycrys talline  ceramic  bodies  could  be 
dramatically  Improved  by  using  powders  having  uniform  particle  sizes.  Also, 
the  time- tempera ture  cycles  needed  for  complete  densif ication  are  lower  in 
temperature  and  shorter  in  time  than  needed  for  conventional  powders.3  The 
superior  microstructural  and  processing  characteristics  result  directly  from 
the  defect-free,  high  coordination  number  green  bodies  that  can  be  produced 
with  uniform  diameter,  agglomerate  free  powders.  Additionally,  the  powders 
must  have  a  small  diameter,  be  equiaxed,  have  high  purity  and  have  specific 
stoichiometry  and  crystallinity. 

Because  existing  powder  synthesis  techniques  could  not  produce  SijNu  and 

SiC  powders  with  these  requisite  characteristics,  we  developed  a  laser  heated 

it  5  6 

gas  phase  synthesis  process.  ’  ’  In  the  laser  process  uniform  nucleation 

and  growth  rates  produce  uniform  powders.  The  program  has  emphasized  the 

synthesis  of  Si,  Si3N4  and  SiC  powders.  The  Si  powders  are  used  for  the 

reaction  bonding  process  In  which  densif ication  and  conversion  to  Si3N4  occur 

simultaneously.  We  have  had  two  specific  objectives.  Our  initial  and 

primary  objective  has  been  to  develop  a  process  for  producing  ceramic  powders 

having  the  desired  characteristics.  To  this  end,  we  have  concentrated  on 

developing  an  analytical  description  of  the  laser  driven  synthesis 
7  8 

process.  *  Secondly,  we  have  begun  to  process  the  powders  into 
9  10  11 

bodies.  This  phase  has  concentrated  on  finding  means  to  disperse  SI 

powders  in  nonaqueous  liquids.  Particle  size  distributions,  agglomerate 
sizes,  zeta  potentials,  and  coagulation  rates  have  been  measured.  Coulombic 
and  van  der  Waals  Interparticle  forces  have  been  calculated. 

Based  on  results  with  metals  and  other  considerations,  we  anticipate 
that  unusual  and  superior  materials  properties  can  be  realized  by  rapidly 
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solidifying  ceramic  materials.  Unusual,  attractive  densifica tion 
characteristics  have  already  been  observed11  with  the  amorphous  powders  that 
can  be  produced  by  the  laser  heated  gas  phase  synthesis  process.  We  have 
developed  the  means  to  produce  rapidly  solidified  ceramics  and  have  begun 
preliminary  characterization  of  these  materials.12  Initially,  these 
experiments  are  focusing  on  materials  that  exhibit  interesting  electrical 
properties.  Also,  prior  to  being  dropped  into  the  quenching  apparatus,  the 
melts  can  be  characterized  to  yield  surface  tensions,  densities  and 
viscosities.  Our  initial  research  has  concentrated  on  an  improved  analytical 
procedure  for  pendant-drops  needed  for  melts  having  surface  tension  and 
densities  characteristic  of  oxides.13 

This  report  consists  of  five  technical  chapters.  Four  are  preprints  of 
papers  that  have  been  submitted  for  publication  and  the  fifth  is  a  statement 
of  the  present  status  of  our  ability  to  disperse  nonionic  powders  in 
nonaqueous  media.  The  principal  contributors  are  reflected  by  the 
authorship. 

Studies  of  the  synthesis  process  have  addressed  many  technical  and 
economic  issues.  The  second  chapter,1'4  entitled  "Growth  of  Precisely 
Controlled  Powders  from  Laser  Heated  Gases",  summarizes  the  general  features 
of  the  synthesis  process  as  it  is  presently  configured,  the  characteristics 
of  resulting  powders  and  projected  manufacturing  costs.  The  results  are 
extremely  encouraging.  The  third  chapter,  entitled  "Powder  Temperature,  Size 
and  Number  Density  Evolution  in  Laser  Driven  Reactions",  reports  the 
procedures  developed  to  monitor  nucleation  and  growth  of  particles  in  laser 
heated  reactants.  This  unique  experimental  technique  provides  a  basis  for 
analysing  reaction  kinetics  as  a  function  of  reactant  partial  pressures  and 
temperatures  with  a  0.5  mm  spacial  resolution.  The  fourth  chapter,  entitled 
"Homogeneous  Nucleation  and  Growth  of  Silicon  Powder  from  Laser  Heated  SiH^”, 
applies  this  diagnostic  technique  to  the  nucleation  and  growth  kinetics  of  Si 
particles  from  SIH^  gas.  Rate  controlling  steps  for  nucleation  and  growth 
steps  were  identified.  The  fifth  chapter,  entitled  "Dispersion  of  Silicon  in 
n-Propanol" ,  summarizes  the  dispersion  results  for  Si  powders.  Based  on 
these  results  and  analyses,  the  feasibility  of  successfully  employing 
coulombic  and  steric  stabilization  are  discussed.  On  the  second  subject,  the 
sixth  chapter,  entitled  "Surface  Tensions  of  Alumina-Containing  Liquids", 
presents  the  first  results  of  characterizing  high  temperature  liquids. 
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Surface  tensions  of  alumina-containing  melts  are  reported  as  a  function  of 
composition  and  ambient  atmosphere. 

The  results  of  characterizations  continue  to  demonstrate  that  the  laser 
heated  gas  phase  powder  synthesis  process  is  capable  of  producing  Si,  Si  ?:4 
and  SiC  powders  having  the  desired  characteristics.  The  particles  can  be 
small,  uniform  in  diameter,  weakly  agglomerated,  spherical,  of  high  pur  ity 
and  controlled  stoichiometry  and  crystallinity.  While  the  maximum  diameter 
achieved  with  Si  (~  2500A)  exceeds  our  target  of  1500A,  to  date  we  have  not 
reached  this  large  a  diameter  with  SijN4  and  SiC.  By  manipulating  process 
variables  we  have  suppressed  nucleation  to  permit  maximum  diameters  to  be 
increased  to  ~  700£;  further  processing  research  is  required  to  achieve 
larger  Si3N4  and  SiC  powders.  Recent  results  with  Si  have  shown  that  hard 
agglomerates  can  form  both  by  sintering  and  melting  mechanisms  when  heated 
excessively.  Thorough  characterization  of  these  types  of  particles  has 
resolved  some  anamolies  between  measured  and  BET  equivalent  diameters  of  the 
larger  diameter  powders.  Overall,  the  process  still  appears  capable  of 
satisfying  the  technical  requirements  and  is  scalable  to  higher  rates  of 
production. 

A  simplified  manufacturing  cost  analysis  was  undertaken  to  provide  a 
first  estimate  of  anticipated  production  costs  and  to  reveal  critical  cost 
factors.  These  results  show  that  the  process  can  produce  submicron  diameter 
powders  at  a  cost  of  $2.60  to  $3.70  per  kg  exclusive  of  reactant  materials. 
The  future  cost  of  SiH4  dominates  anticipated  manufacturing  costs.  Even  with 
higher  than  anticipated  costs  for  SiH4  20$/kg),  this  process  appears 
capable  of  producing  submicron  SiC  powders  at  substantially  lower  costs  than 
those  produced  by  grinding  Acheson  material.  This  very  important  result  is 
the  basis  for  scaling  the  process  to  pilot  scale  production  levels. 

The  direct  measurement  of  nucleation  and  growth  kinetics  formed  an 
important  link  between  our  previously  developed  empirical  relationships 
between  process  variables  and  powder  characteristics  and  the  existing  kinetic 
models  for  low  growth  rate  CVD  processes.  The  new  diagnostic  technique  is 
based  on  simultaneous  HeNe  laser  light  scattering  and  transmission 
measurements.  The  results  permit  calculation  of  the  local  emissivity,  number 
density  and  particle  diameter  in  the  reaction  zone. 

Combined  with  measurements  of  brightness  temperatures  and  calculated  gas 
velocities,  nucleation  and  growth  rates  are  derived  as  a  function  of 
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temperature  and  reactant  partial  pressures.  The  calculated  results  are 
extremely  sensitive  to  the  assumed  complex  index  of  refraction.  While  this 
introduces  a  potential  for  error,  the  sensitivity  can  be  used  to  interpret 
reaction  mechanisms  when  calculated  results  are  compared  with  the  results  of 
other  independent  characterizations. 

The  diagnostic  technique  gives  heating  rates  (106-10e°C/s) ,  maximum 
temperature,  dwell  time  at  maximum  temperature  (~  10~4s),  and  cooling  rates 
(~  10-5°C/s).  For  Si  powder,  it  also  provided  direct  evidence  that  the 
nucleation  process  stops  early  in  the  reaction  (~  27.  conversion),  that 
agglomerates  do  not  form  in  the  heated  zone  and  that  the  Si  particles  are 
amorphous  when  formed  and  then  may  crystallize  if  heated  sufficiently.  The 
growth  processes  could  be  monitored  from  a  minimium  particle  diameter  of 
~  80*. 

The  nucleation  and  growth  kinetics  of  Si  particles  from  SiH^  gas  were 
studied.  Particles  nucleate  homogeneously  at  rates  in  excess  1014  per 
cm3/s.  The  high  level  of  initial  supersaturation  presents  no  barrier  to 
nucleation,  rather  the  nucleation  rate  appears  controlled  by  the 
concentration  of  Si  atoms  in  the  vapor  phase.  Growth  rate  is  controlled  by 
vapor  phase  transport  of  reactants  to  the  particle  surfaces.  The 
condensation  coefficient  is  estimated  to  be  0.002-0.02.  Particle  growth 
terminates  when  reactants  are  depleted;  the  final  particle  size  depends  on 
the  initial  nuclei  density  and  the  reactant  concentration.  The  separation  of 
the  nucleation  and  growth  processes  intrinsically  tends  to  yield  uniform 
powders. 
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Building  on  earlier  research,  the  dispersibility  of  Si  powders  in 

n-propanol  has  been  studied.  This  work  incorporated  a  more  complete 
characterization  of  the  constituent  powders.  Particle  size  distribution, 
agglomeration,  agglomerate  size,  bulk  and  surface  purities  and  (-potential 
were  characterized  for  pure  and  boron  doped  powders.  Dispersion  stability 
and  agglomeration  rates  were  evaluated  as  functions  of  solids  and  water 
contents.  The  results  were  interpreted  in  terms  of  DLVO  theory.  Stable 
dispersions  can  be  produced  with  low  solids  content  and  dry  n-propanol.  The 
low,  dispersed  coulombic  potential  barrier  resulting  from  a  large  double 
layer  does  not  provide  adequate  repulsion  for  close  lnterpartlcle  spacing. 
Stabilization  by  means  of  ionic  additives  which  compress  the  double  layer  may 
be  possible.  Alternatively,  steric  barriers  ccn  be  used  to  counteract  van 


der  Waals  attractive  forces  at  close  interparticle  spacings.  Because  the 
theory  for  nonaqueous  dispersions  is  not  well  developed,  the  dispersion 
research  will  of  necessity  proceed  as  a  largely  empirical  study  guided  by 
theory. 

The  laser  heated  floating  zone  crystal  growth  apparatus  was  modified  to 
permit  melt  surface  tension  and  rapid  solidification  experiments.  Both 
static  and  dynamic  techniques  are  to  be  used  to  provide  independent 
measurements  of  surface  tension  to  density  ratios.  The  dynamic  technique 
will  also  permit  measurement  of  melt  viscosity.  A  piston  and  anvil  quenching 
apparatus  has  been  installed  to  splat  quench  high  temperature  liquids.  Melts 
from  Y203~Zr02,  Al203-Zr02>  Al203~Ti02,  Al203~Cr203,  Al203~Mg0,  Sc203~Ta205 
and  Li20-B203  systems  have  been  quenched  and  characterized. 

The  results  of  the  static,  pendant-drop  surface  tension  measurements  are 
reported  in  detail.  An  improved  analytical  procedure  based  on  comparing 
calculated  and  observed  melt  contours  was  developed.  This  new  technique 
permits  acceptable  experimental  accuracy  with  the  short  maximum  pendant-drop 
lengths  achievable  with  the  oxide  melts.  Surface  tensions  of  A1„C>3 
containing  melts  with  Cr203 ,  MgO,  Ti02  and  Zr02  additives  were  measured  in 
air,  helium  and  helium-hydrogen  atmospheres.  The  surface  tensions  tended  to 
increase  with  oxygen  partial  pressure.  In  air,  the  effect  of  additives  on 
surface  tension  parallelled  their  effect  on  the  melting  point.  In  most 
cases,  the  compositional  effects  were  approximately  linear;  TI02,  an 
exception,  caused  a  sharp  decrease  in  y  for  small  additions  that  were 
followed  by  essentially  constant  y's  for  further  additions.  The  effect  of 
additives  in  neutral  and  reducing  atmospheres  was  more  complex  and  did  not 
follow  specific  trends. 

During  the  first  year,  the  research  has  addressed  all  topics  in  the 
contractual  scope  of  work.  Technically  and  economically,  the  laser  heated 
gas  phase  synthesis  process  appears  to  be  a  viable  means  of  producing  Si, 
SijN^  and  SIC  powders.  Many  Issues  remain  to  be  resolved  before  the  highly 
perfect  microstructures  achieved  with  TiC>2  powders  are  realized  with  these 
nonionic  materials;  but,  they  appear  tractable.  The  rapidly  solidified 
ceramics  continue  to  be  interesting  based  on  results  of  microstructural 
characteristics. 
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CHAPTER  II 


GROWTH  OP  PRECISELY  CONTROLLED  POWDERS  PROM  LASER  HEATED  GASES 

by 

J.  S.  Haggerty 

ABSTRACT 

Powders  having  superior  characteristics  have  been  synthesized  from 
laser  heated  gases.  These  superior  properties  result  from  the  unusually 
well  controlled  process  conditions  that  are  achievable  with  laser  heating. 
Manufacturing  costs  appear  lower  than  conventional  submicron  powders. 

II. 1.  INTRODUCTION 

This  paper  will  summarize  the  principal  features  of  a  powder  synthesis 
process  developed  at  M.I.T.  Many  of  the  technical  details  have  already  been 
described1-6  and  others  are  being  published.  The  latter  include  papers 
based  on  the  theses  of  Flint7  and  Marra8,  as  well  as  contributions  by 
Suyama9,  Golino10  and  Gattuso11.  This  synthesis  process  was  developed 
specifically  to  enable  the  making  of  powders  having  characteristics 
considered  ideal  for  processing  into  ceramic  parts  has  been  discussed  by 
Bowen.12  These  ideal  powder  characteristics  include:  uniform  particle 
size,  small  particle  size,  equiaxed  shapes,  freedom  from  agglomerates,  and 
absence  of  impurities  or  deviations  from  desired  stoichiometry. 

Our  processing  research  has  concentrated  on  the  nonoxide  particulate 
materials  Si,  Si3N^  and  SiC;  to  a  lesser  extent  oxides  have  also  been  made 
by  this  process.  Based  on  our  results,  it  is  apparent  that  laser  heated  gas 
phase  processes  can  be  employed  with  many  different  materials  to  make  either 
powders  or  monolithic  coatings. 

A  gas  phase  process  was  selected  for  making  powders  having  ideal 
characteristics.  Experience  with  Ti02  powder  synthesis  shows  that  gas  phase 
processes  can  produce  inexpensively  powders  having  the  desired  characteris¬ 
tics.  Heating  techniques  that  were  considered  included  rf  and  microwave 
plasmas,  arcs,  heated  tubes,  reactive  sputtering  techniques  and  optical 
heating  means.  The  latter  were  elected  because  they  appeared  to  combine 
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many  desirable  features.  Although  broadband  optical  sources  can  be  used  to 
heat  gases,  they  are  typically  inefficient  because  the  gases  have  few, 
narrow  absorption  lines.  Heating  with  a  laser  can  be  very  efficient  with 
some  lasers  when  emission  and  absorption  lines  coincide.  The  general 
process  characteristics  achievable  with  laser  heating  that  were  considered 
important  include:  a  small,  well  defined  heating  zone  can  be  created,  all 
gas  molecules  can  be  subjected  to  similar  time- temperature  histories, 
process  variables  can  be  manipulated  in  reproducible  and  defined  manners, 
controlled  atmosphere,  cold-wall  reaction  vessels  are  straight  forward. 

Specifically,  we  have  used  C02  lasers  as  heat  sources  for  the  gas  phase 
reactions.  This  type  of  laser  emits  on  approximately  100  different 
wavelengths  between  nominally  9.2  and  11.0  pm.  Several  of  these  lines  have 
been  investigated  in  our  research;  fortunately,  many  of  the  useful  reactant 
gases  exhibit  strong  absorp tivities  for  the  P(20)  line  at  10.591  pm.  This 
is  the  highest  gain  line  exhibited  by  a  CO^  laser  and  is  both  the  most 
efficient  choice  and  is  the  line  usually  emitted  by  the  simple  untuned 
lasers.  These  lasers  convert  electricity  to  light  with  approximately  15!. 
efficiency;  better  than  many  alternative  heating  means.  This  laser 
equipment  is  extremely  reliable,  controllable  and  has  a  reasonable  cost. 

Reactant  gases  are  heated  by  absorbing  IR  photons  in  these  laser 
processes.  They  differ  from  photochemical  reactions  induced  by  ultraviolet 
or  visible  wavelength  photons  because  absorption  of  single  or  multiple  IR 
photons  generally  leaves  the  molecules  in  vibrationally  or  rotationally 
excited  electronic  ground  states  where  they  ultimately  dissociate  only  when 
sufficient  energy  has  been  absorbed.  With  the  high  gas  pressures  used, 
these  reactions  basically  proceed  as  thermal  reactions  with  a  Boltzmann 
distribution  of  energies.  Even  if  collisions  and  other  relaxation  processes 
occur,  these  reactions  may  differ  from  conventional  reactions  because  the 
short  times  and  rapid  heating  rates  may  result  in  different  intermediate 
states.  Our  investigation  has  not  presumed  unusual  reaction  paths  and  has 
proceeded  on  the  premise  that  these  laser  induced  reactions  will  differ  from 
conventional  reactions  in  the  level  of  process  control  and  the  process 
conditions  that  can  be  achieved  uniquely  with  them  and  that  they  will  have 
value  because  of  the  superior  properties  of  the  resulting  powders. 
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II. 2.  LASER  HEATED  POWDER  SYNTHESIS  PROCESS 


II. 2. A.  Process  Description 

Powders  have  been  synthesized  with  crossflow  and  counterflow  gas 
stream  -  laser  beam  configurations  and  with  both  static  and  flowing  gases. 
Most  of  the  process  research  has  been  conducted  with  the  reaction  cell  shown 
schematically  in  Figure  II. 1. 

In  a  crossflow  configuration,  the  laser  beam  having  a  Gaussian  shaped 
intensity  profile  orthogonally  intersects  the  reactant  gas  stream  possessing 
a  parabolic  velocity  profile.  The  laser  beam  enters  the  cell  through  a  KC1 
window.  The  premixed  reactant  gases,  under  some  conditions  diluted  with  an 
inert  gas,  enter  through  a  1.5  mm  ID  stainless  steel  nozzle  located  2-3  mm 
below  the  laser  beam.  A  coaxial  stream  of  argon  is  used  to  entrain  the 
particles  in  the  gas  stream.  Argon  is  also  passed  across  the  inlet  KC1 
window  to  prevent  powder  build-up  and  possible  window  breakage.  Cell 
pressures  are  maintained  between  0.08  to  0.90  atm  with  a  mechanical  pump  and 
throttling  valve.  At  present,  the  powder  is  captured  in  a  microfiber  filter 
located  between  the  reaction  cell  and  vacuum  pump.  In  the  future,  the 
particles  will  be  collected  by  electrostatic  precipitation  or  by  separation 
with  a  fluid  used  for  dispersion. 

To  insure  uniform  laser  intensities  throughout  the  reactants,  the  gas 
stream  has  been  made  optically  thin.  Typically  only  2-15%  of  the  incident 
light  is  absorbed  by  the  reactants  with  this  experimental  apparatus 
configuration.  Obviously  more  efficient  use  of  the  expensive  laser  energy 
is  essential  for  a  production  process.  The  counterflow  configuration 
permits  the  simultaneous,  complete  capturing  of  the  light  energy  and  uniform 
time  -  temperature  histories.  Although  it  is  one  alternative  that  probably 
satlslfies  production  requirements,  the  counterflow  apparatus  is  not  well 
suited  to  an  experimental  program  in  which  process  variables  are  manipulated 
over  orders  of  magnitude.  There  are  several  other  equipment  configurations 
that  capture  all  of  the  light  while  achieving  uniform  heating.  These 
Include  multiple  radially  opposed  beams,  reflecting  the  beam  repeatedly 
through  an  optically  thin  gas  stream  as  we  are  doing  in  the  laser  heated 
chemical  vapor  deposition  process13,  and  intersecting  an  elongated,  strip 
shaped  gas  stream  at  an  oblique  angle  of  Incidence  with  one  or  more  laser 
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l.  Powder  synthesis  cell. 3 


TABLE  II. 1 

PROCESS  CONDITIONS  USED  FOR  SI,  SljN^  AND 


Process  Variable 


Si  Si3N4 


Cell  Pressure  (atm) 

SiH4  Flow  Rate  (cc/min) 
NH3  Flow  Rate  (cc/min) 
C2H4  Flow  Rate  (cc/min) 
Ar  Flow  Rate  (Annulus 
plus  window)  (cc/min) 
Laser  Intensity  (W/cm2) 
Reaction  Zone  Temp.  (°C) 


0.2-0.  9  0.2-0.  9 

5.4-110  5.4-40 

0  44-110 

0  0 


560-1450  1000-1100 
176-5. 4xl03  530-lxl05 
750-1390  675-1390 


SIC 


SiC 


0.2-0.  8 
11-45 
0 

9.0-45 

1000 

530-5.2x103 

865-1930 


beams.  The  virtually  complete  conversion  of  laser  light  to  process  heat 
appears  quite  feasible. 

Principally,  Si,  Si^N^  and  SiC  powders  have  been  made  from  appropriate 
combinations  of  SiH^,  NH^,  and  C2H^  gases.  CH^  was  investigated  but  we 
found  that  it  did  not  react  to  form  SiC  with  the  laser  heated  process 
conditions.  Alternative  reactant  gases  such  as  CH3*SiH3  and  Cl2SiH2  have 
been  investigated  for  SiC;  others  have  also  been  used  in  preliminary 
experiments  to  make  Al^Oj,  TiC>2  and  B^C.  Diborane  (B2H6)  was  used  to  boron 
dope  Si  and  SiC  powders  acting  as  a  sintering  aid. 14  Thus,  while  most  of 
this  process  research  has  focused  on  a  limited  set  of  compounds  and 
reactants,  it  is  apparent  that  laser  induced  reactions  are  applicable  to  a 
broad  range  of  materials. 

The  specific  process  conditions  employed  to  make  the  powders  have  been 
manipulated  over  an  extensive  range.  The  investigated  conditions  are 
summarized  in  Table  II. 1.  More  recent  synthesis  runs  have  employed  higher 
pressures,  higher  mass  flow  rates  and  near  zero  argon  flow  rates  for  the 
window  stream.  Fortunately,  these  latter  conditions  combine  more  desirable 
powder  characteristics  with  economically  more  attractive  process 
parame ters. 

Heating  rates,  nucleation  rates,  growth  rates  and  temperatures  can  be 
estimated  surprisingly  accurately  based  on  measured  zone  and  particle 
dimensions,  calculated  gas  velocities  and  uncorrected  pyrome trically 
determined  temperatures.  For  typical  reaction  conditions  employing  a  total 
reactant  gas  flow  rate  of  approximately  100  cm3/min  and  a  pressure  of 
0.2  atm,  the  velocity  of  the  gas  decreases  from  approximately  500  cm/s  at 
the  nozzle  to  350  cm/s  at  the  center  of  the  laser  beam.  With  a  NH^/SiH^ 
flame,  the  reaction  commences  approximately  3  to  5  mm  into  the  laser  beam. 
Thus,  the  exposure  time  needed  to  initiate  the  reaction  is  nominally 
10~3  s.  For  a  reaction  temperature  of  approximately  1000°C,  this  indicates 
a  gas  heating  rate  of  approximately  106#C/s.  Time- temperature  profiles  have 
been  computed3  as  a  function  of  varied  process  parameters  using  calculated 
velocities,  actual  laser  intensities  and  absorptlvltles  that  we  have 
measured.  Qualitatively,  the  calculations  agree  with  observed  effects. 

For  most  process  conditions,  the  reactants  are  essentially  completely 
consumed;  this  is  particularly  important  with  expensive  reactants.  Also, 
particle  growth  is  limited  by  depletion  of  the  reactants  rather  than  elapsed 
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TABLE  II. 2 


l 

SUMMARY  OF  POWDER  CHARACTERISTICS 


Powder  Type 

SI 

Si3N4 

SIC 

Mean  Diameters  (') 

Standard  Deviation  of  Diameters 

190-1700 

46 

75-500 

23 

200-500 

-25 

(%  of  Mean) 

Impurities  02(wt.  %) 

0.06-  0.7 

-0. 3 

0.3  -  1.3 

Total  Others  (ppm) 

<200 

<200 

h'A 

Major  Elements 

Ca , Cu , Fe 

Al.Ca 

NA 

Stoichiometry 

— 

0-60t 

o-i  or. 

(excess  SI) (excess  C  or  SI) 

Crystalllnl ty 

crys  tal 1 Ine 

amorphous 

crys  talllne 

-amorphous 

-crystalline 

SI  and  SIC 

Grain  Slze:Mean  Diameter 

1/5-1/3 

-1/2 

1/2-1 

14 


growth  time  which  permits  particle  size  to  be  manipulated  through  control  of 
the  nucleation  process. 


II. 2. B.  Powder  Characteristics 

The  Si,  Si3N4  and  SiC  powders  all  exhibit  the  same  general  features  and 
appear  to  match  the  idealized  characteristics  we  sought.  Table  II. 2 

summarizes  the  range  of  properties  that  have  been  achieved  with  these 

3  9 

particulate  materials.  ’ 

A  bright  field  TEM  photomicrograph  of  a  typical  SiC  powder  is  shown  in 
Figure  II. 2.  The  particles  are  spherical  and  uniform  in  size.  The  SisN 
and  SiC  powders  are  usually  smaller  and  have  a  narrower  size  distribution 
than  the  Si  powders.  This  micrograph  also  indicates  the  powders  were 
observed  in  chainlike  agglomerates  after  being  captured  in  a  filter.  Neck 
formation  has  been  observed  by  TEM  between  Si  particles,  but  not  with  the 
Si3N’u  and  SiC  powders.  Dispersion  results15  have  shown  primary  bonding  does 
not  exist  in  necks  between  Si  particles  because  light  scattering  and  photon 
correlation  spectrometer  characterizations  of  the  sols  showed  that  the 

dispersed  particle  sizes  are  equal  to  those  of  the  individual 
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particles.  ’  ’  It  is  presumed  that  the  same  will  be  true  of  Si;S'4  and 

SiC  powders  that  do  not  exhibit  interparticle  necks. 

The  BET  equivalent  spherical  diameter  and  the  diameter  measured  from 
TEM  micrographs  have  always  been  nearly  equal.  This  indicates  that  the 
particles  have  smooth  surfaces,  no  porosity  accessible  to  the  surface, 
relatively  narrow  size  distributions,  and  nearly  spherical  shapes.  Powder 
densities,  measured  by  He  pycnometry,  indicated  the  particles  had  no 
Internal  porosity. 

Chemical  analyses  indicate  that  the  oxygen  content  is  generally  less 
than  1.0%  by  weight  and  some  powders  are  as  low  as  0.05  wt.%.  The  total 
cation  impurities  are  typically  less  than  200  ppm.  For  SijN4  and  SiC 
powders,  the  stoichiometry  varied  substantially  depending  on  the  process 
conditions.  For  Si3N4  powder,  the  stoichiometry  ranged  from  nearly  pure 
SijN4  (<  1,0  wt.%  excess  Si)  to  SijN4  +  60  wt.%  excess  Si.  Near 
stoichiometric  values  are  caused  by  increased  laser  intensity,  increased 
pressure,  and  lower  gas  stream  velocities.  C2H4/SiH4  runs  produced  SiC 
powders  with  up  to  10  wt.%  excess  C  or  Si.  The  more  stoichiometric  powders 
were  produced  with  increased  laser  intensity  and  close  to  s toichiometric 


reactant  gas  ratios.  In  both  the  Si^N^  and  SiC  powders,  it  appears  the 
excess  Si  and  C  are  distributed  uniformly  throughout  the  individual 
particles. 

Except  when  synthesized  under  very  low  laser  intensities,  the  Si 
powders  are  crystalline.  In  all  cases,  the  crystallite  size  was  a  fraction 
(1/5  -  1/3)  of  the  BET  equivalent  diameter,  indicating  the  individual 
particles  are  polycrystalline.  Marra  found8  that  the  particle  size  to  grain 
size  ratio  reflected  the  nucleation  and  growth  of  crystals  in  the  amorphous 
particles.  Virtually  all  process  conditions  produced  Si^N^  powder  which 
was  amorphous.  High  laser  intensities  and  high  reactant  pressure  resulted 
in  stoichiometric,  crystalline  SijN^  powders.  All  SiC  powders  were 
crystalline.  Most  process  parameters  produced  polycrystalline  SiC  particles 
with  the  BET  equivalent  spherical  diameter  approximately  twice  as  large  as 
the  crystallite  size  measured  by  X-ray  line  broadening. 

Powders  of  other  materials  have  not  been  characterized  as  extensively 
as  these  three  nor  have  the  interrelationships  between  process  variables  and 
powder  characteristics  been  studied  yet.  They  all  have  poorly  developed 
crystalline  structures  or  extremely  small  grain  sizes  based  on  X-ray 
diffraction  analyses  and  have  small  diameters  based  on  TEM  and  BET 
characteristics. 

II. 2. C.  Process  Attributes 

Qualitatively,  it  is  now  clear  that  the  laser  heated  powder  synthesis 
process  works;  in  fact,  it  works  very  well.  It  is  stable,  simple, 
dlagnosable,  efficient  and  the  resulting  powder  characteristics  satisfy  a 
stringent  set  of  requirements.  It  is  the  inherent  stability  and  the  ability 
to  make  real  time  diagnostics  that  permit  high  quality  powders  to  be 
produced,  Insuring  they  have  a  high  value.  The  simplicity  and  both  material 
and  energy  efficiencies  should  permit  costs  to  reach  viable  levels. 

Newly  developed  diagnostics  based  on  light  scattering  and  transmittance 
measurements  permit  nucleation  and  growth  processes  to  be  analyzed  in  terms 
of  temperature  and  reactant  partial  pressures  with  <  1  mm  spacial  resolution 
throughout  the  reaction  zone  and  particulate  plume.  They  also  permit 
important  observations  about  the  crystalline  state  of  the  reaction  product 
and  the  agglomeration  process.  Using  an  experimental  technique  developed  by 
Flint7  and  Marra8  and  modifications7  of  computational  techniques  reported  by 


Sarofim^7,  it  is  possible  to  calculate  the  local  emissivity,  particle 

diameter  and  number  density  from  the  extinction  and  scatter-extinction  ratio 

results.  Mass  balance  calculations,  based  on  local  number  density  and 

particle  diameter,  permit  the  local  reactant  concentrations  to  be 

calculated.  These,  combined  with  local  temperatures,  permit  reaction 

kinetics  to  be  interpreted  in  terms  of  rate  controlling  models.  They  also 

7  e 

provide  a  bas*s  for  real  time  process  control.  Typical  results  ’  are  given 
in  Figures  II. 3,  II. 4  and  II. 5  for  Si  powders  made  from  SIHj..  The 
temperatures  shown  in  Figure  II. 3  are  calculated  from  the  pyrometrically 
determined  brightness  temperatures  and  the  calculated  emissivi ties . 

These  calculations  employ  the  complex  index  of  the  particulate 
materials  which  are  strongly  dependent  on  crystallinity,  temperature, 
stoichiometry  and  wavelength.  This  variance  can  be  used  to  gain  insights 
about  the  reaction  process.  For  either  amorphous  or  polycrystalline 
indices,  Figure  II. 5,  the  calculated  number  densities  decrease  rapidly  over 
a  distance  of  approximately  1  mm  (~  3  x  10_4s)  and  then  remain  constant. 
Supported  by  direct  STEM  observations  and  annealing  studies®,  these  results 
show  that  the  Si  particles  form  with  an  amorphous  structure  then  crystallize 
when  the  gases  become  hotter  as  they  penetrate  into  the  laser  beam."’  The 
actual  number  density  follows  the  dotted  line  from  the  amorphous  curve 
(~  3  mm)  to  the  polycrystalline  curve  (~  4.5  mm)  and  remains  essentially 
constant  (N  =  4  x  10’i2cm-3)  throughout  the  growth  process.  The  cessation  of 
nucleation  is  an  essential  requirement  for  achieving  the  desired  uniform 
particle  size.  Equally  Important,  this  result  shows  that  agglomeration  does 
not  occur  within  the  hot  region  where  hard  agglomerates  could  form  in 
agreement  with  the  constant  particle  diameter  regions  shown  in  Figure  II. 4. 
This  result  shows  that  the  chainlike  agglomerates  probably  result  from  being 
captured  in  the  filter,  extraction  or  preparation  of  TEM  samples. 

In  close  agreement  with  the  simplified  analyses,  these  direct 
observations  show  that  the  process  is  characterized  by  rapid  heating  rates 
(~  106°C/s),  short  reaction  times  (~  5  x  10-i,s)  and  rapid  cooling  rates 
(~  105°C/s).  The  high  heating  rates  cause  high  nucleation  rates  and 
consequently  small  particle  sizes.  The  short  time  at  high  temperatures 
minimizes  problems  with  forming  hard  agglomerates  when  particles  collide 
through  Brownian  motion.  The  cooling  rates  are  fast  enough  to  quench  non¬ 
equilibrium  structures.  The  maximum  exposure  temperature  can  also  be 
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II. 3.  Temperature  versus  distance  from  nozzle  for  several  Si  synthesis 
runs. ' J 


II. A.  Particle  size  versus  distance  from  nozzle  for  several  synthesis 
runs . 7  ,8 
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II. 5.  Particle  number  density  versus  distance  from  nozzle  for  a  typical  Si 
synthesis  run.7,6 


precisely  controlled  largely  independent  of  heating  and  cooling  rates.  This 
unusual  feature  permits  the  powder's  crystalline  state  to  be  manipulated. 

In  combination,  the  attributes  of  this  process  constitute  a  very  precise 
means  of  producing  ceramic  powders. 

The  costs  of  manufacturing  powders  by  the  laser  heated  processes  have 
been  estimated  and  compared  to  conventional  powders.  The  results  are  very 
encouraging. 

Our  manufacturing  cost  estimates  are  based  on  the  simplified  technique 
described  by  Aries  and  Newton.18  Material,  utility,  energy  and  labor  costs 
are  treated  conventionally.  The  major  simplification  is  that  the  plant  cost 
is  assumed  to  be  a  multiple  of  the  equipment  cost  and  that  annualized  plant 
cost  is  0.202  x  total  plant  cost.  Chemical  engineering  experience  shows 
that  the  "plant  cost”  multiple  should  be  4.74  for  a  fluid  process;  lower 
multiples  are  used  for  simpler,  non  hermetic  processes.  Labor  costs  are 
treated  as  1.79  x  direct  salaries  and  annual  building  costs  as  0.052  x  total 
building  costs. 

This  analysis  was  based  on  a  single  gas  tip  since  this  would  be  the 
scaling  unit  for  a  plant.  It  presumed  most  of  the  same  operating  conditions 
we  employ  with  the  laboratory  scale  process;  they  are  summarized  in 
Table  II. 3.  The  assumed  2.0  cm  diameter  gas  tip  is  the  only  significant 
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TABLE  II. 3 


MANUFACTURING  RATES  AND  COSTS 

•  Assumptions 

Gas  tip  diameter  =  2.0  cm 

Gas  velocity  =  500  cm/sec 

Gas  pressure  =  1  atm 

Reaction  efficiency  =  100  % 

Capture  efficiency  =  100% 

Operational  efficiency  =  91%  (8000  hr/yr) 

Heating  - 

0  -  700°C  by  resistance 
700  -  1000°C  by  laser 
no  heat  recovery 


Resul ts 

Si 

SiC 

Production  rate:  kg/hr 

7.2 

5.1 

6.8 

ton/yr 

57.6 

40. 8 

5A.4 

Laser  power  requirement  (kW) 

1.5 

1.3 

1.5 

Total  process  energy  (kWhr/kg) 

1.8 

2.1 

1.9 

Costs  ($/kg) 


Labor 

0.77 

1.09 

0.82 

Plant 

1.66 

2.34 

1.76 

Building 

<0.01 

<0.01 

<0.01 

Utilities 

0.19 

0.25 

0.21 

Raw  materials 

2.28-22.80 

1.62-13.95 

1.77-16. 17 

Total 

A. 29-25. 44 

5.32-17.65 

4.57-18.97 

2 


departure  from  current  practice.  This  value  is  approximately  the  maximum 

laser  beam  penetration  depth  that  can  produce  uniform  local  intensities  with 

multiple,  radially  apposed  laser  beams.  The  assumed  reaction  and  powder 

capture  efficiencies  reflect  our  current  values.  We  presumed  that  preheated 

gases  would  be  injected  into  the  laser  beam  to  minimize  unnecessary  laser 

heating;  but,  no  heat  recovery  was  assumed.  The  plant  was  scheduled  for 

full  time  operation  (8760  hrs/yr  maximum)  and  assumed  to  actually  produce 

91%  of  available  time  (8000  hrs/yr).  Four  full  shifts  having  one  person  for 

three  tips  was  assumed.  Current  bulk  costs  for  NH?  (0.52  $4/kg)  srd 

(0.5  $/kg)  were  used  for  these  materials.  Silane  costs  were  estimated 
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informally  from  several  sources.  *  ’  The  cost  analysis  was  done  based 

on  a  silane  cost  of  2.0  and  20.0  $/kg.  A  figure  of  10-15  $/kg  appears 
reasonable  for  silane  with  somewhat  higher  impurity  levels  than  semi¬ 
conductor  grades. 

The  results  of  the  manufacturing  cost  analysis  are  given  in  Table  II. 3. 
They  show  that  the  assumed  process  configuration  is  well  matched  to 
available  C02  laser  equipment.  The  presumed  gas  tip  diameter  and  gas  flow 
rate  require  lasers  in  the  1.5  kW  range.  Also,  the  output  per  tip  is  40-57 
tons  per  year;  a  level  that  is  well  into  a  viable  manufacturing  scale. 
Finally,  calculated  manufacturing  costs  range  from  4.30  to  25.00  $/kg 
depending  mostly  on  the  silane  cost.  The  dominance  by  raw  material  costs 
will  be  ‘rue  of  other  gas  phase  processes  using  these  materials.  Their 
distinctions  will  be  the  quality  of  the  resulting  powders.  These  costs  pose 

questions  regarding  the  accuracy  of  the  analysis,  the  opportunities  for  cost 

reduction  and  the  relative  costs  of  conventional  materials. 

Generally  we  feel  that  the  cost  projection  is  conservative.  The  only 
factor  not  considered  is  maintenance.  The  "plant  cost"  estimating  factor 
appears  too  high  because  the  required  equipment  cost  is  made  up  largely  of  a 
single,  easily  installed  laser.  For  instance,  equipment  having  the  optics 
and  laser  power  assumed  in  this  analysis  was  just  designed  and  constructed 
at  M.I.T.  on  a  first  time  basis  for  much  less  than  half  of  the  ~  $470,000 
figure  used  in  the  calculation.  We  believe  this  cost  will  be  less  than 

indicated.  More  importantly,  there  are  many  opportunities  for  reducing  the 

raw  material  costs  made  possible  by  changing  chemistries.  Our  initial 
choice  of  silane  as  a  silicon  source  made  no  consideration  of  cost.  A  major 
task  In  our  new  "scale  up"  research  program  will  be  the  consideration  of 
alternative  reactants. 
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II. 6  Sales  price  in  5000  kg  lots  for  Acheson  SIC  powder22,  specific 

comminution  energy23  and  manufacturing  cost  for  SIC  powders  made 
from  laser  heated  gases  as  functions  of  particle  size. 


Figure  II. 6  summarizes  present  sales  prices  of  Acheson  SiC  in  5000  kg 
lots  as  a  function  of  particle  size.22  The  price  is  essentially  constant 
for  powder  sizes  greater  than  100  urn.  Below  100  pm,  the  cost  increases 
rapidly  with  decreasing  particle  size.  A  1/d  extrapolation  was  used  to 
project  prices  into  the  submicron  range  needed  for  sintering.  Also  shown  In 
Figure  II. 6,  Is  a  generalized  curve  of  comminution  energy  per  unit  mass  as  a 
function  of  particle  size.23  For  coarse  powders,  it  Is  independent  of 
particle  size;  for  fine  and  ultrafine  powders  it  usually  increases  with 
decreasing  particle  size  at  a  rate  between  1/d  and  1/d2.  The  similarity 
between  the  two  curves  Is  striking;  although,  it  should  be  noted  that  the 
values  of  these  specific  comminution  energies  constitute  a  small  fraction  of 
the  sales  prices.  The  manufacturing  costs  for  powders  made  by  a  comminution 
process  must  follow  curves  having  shapes  similar  to  those  shown  in 
Figure  II. 6  and  also  they  cannot  be  displaced  very  much  for  profit  factors. 
The  "coarse",  constant  cost  is  largely  set  by  energy  and  handling  costs,  the 
diameter  dependence  usually  begins  between  5  and  100  ym,  below  that,  costs 
increase  at  rates  between  1/d  and  1/d2.  These  powders  will  have  the  typical 
broad  particle  size  distribution  and  contaminations  resulting  from 
comminution  processes. 


A. 


Our  SIC  manufacturing  cost  estimate  range  is  shown  in  Figure  II. 6  for 
comparison  with  commercial  powders.  These  costs  are  independent  of  particle 
diameter  because  controlling  the  nucleation  process  has  no  obvious 
consequence  to  equipment  costs  or  mass  flow  rates.  Even  with  ;  robably 
higher  than  actual  cost  figures,  the  laser  heated  gas  phase  synthesis 
process  has  a  much  lower  manufacturing  cost  than  the  sales  price  of  present 
submicron  SIC  powders.  Also,  these  powders  have  a  narrow  size  distribution 
and  are  free  of  unwanted  contaminates. 

11. 3.  CONCLUSIONS 

The  conclusions  are  straight  forward  and  are  very  optimistic.  A  new 
powder  process  has  been  developed.  It  is  simple,  stable,  controllable  on  a 
real-time  basis  and  applicable  to  many  materials.  The  characteristics  of 
the  resulting  powders  satisfy  criteria  for  an  ideal  powder  and  individual 
properties  can  be  manipulated.  Projected  manufacturing  costs  are  lower  than 
conventional  powders.  It  appears  to  us  that  this  new  process  is  capable  of 
producing  better  materials  at  lower  cost  than  are  presently  available. 
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CHAPTER  III 


POWDER  TEMPERATURE,  SIZE,  AND  NUMBER  DENSITY 
IN  LASER  DRIVEN  REACTIONS 
by 

J.  H.  Flint,  R.  A.  Marra,  and  J.  S.  Haggerty 

ABSTRACT 

A  technique  to  accurately  measure  the  temperature  of  powders  in  a  laser 
driven  reaction  has  been  developed.  Particles  are  formed  by  heating 
reactant  gases  with  a  150  watt  C02  laser.  The  brightness  temperature  of  the 
particulate  cloud  was  measured  with  a  micro-optical  pyrometer.  The 
emissivity  was  determined  from  scattering  and  transmission  measurements.  A 
correction  for  high  turbidity  was  derived.  The  scattering  and  transmission 
measurements  also  allow  the  determination  of  the  size  and  number  density  of 
the  particles.  The  temperature  and  particle  size  as  a  function  of  height 
are  reported  for  five  silicon  powder  reaction  conditions,  and  for  one 
silicon  carbide  reaction. 

The  measurements  indicate  that  particles  are  often  nucleated  before  the 
reactant  gas  has  reached  the  C02  laser  beam.  The  reaction  zone  temperature 
decreases  once  most  of  the  reactant  gas  is  consumed,  unless  the  produced 
powder  absorbs  10.6  um  radiation  as  does  SiC.  The  silicon  particles 
nucleate  as  amorphous  silicon,  and  then  crystallize  as  they  move  into  hotter 
regions  of  the  reaction  zone. 

III.l.  INTRODUCTION 

A  method  of  producing  sinterable  powders  through  the  heating  of  gases 

with  a  C02  laser  has  been  developed  which  has  potential  application  in  the 

manufacture  of  silicon  nitride  and  silicon  carbide  as  well  as  other  high 

12 

temperature  ceramics.  Although  variation  in  the  product  powder 
characteristics  have  been  studied  extensively  as  a  function  of  process 
parameters,  the  properties  of  the  reaction  itself  have  been  more  difficult 
to  characterize.  The  most  important  of  these  properties  are:  the  thermal 
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range  within  which  the  reaction  takes  place;  and  the  evolution  in  tine  (and 
space)  of  the  particle  size,  particle  number  density,  and  crystallinity. 

This  paper  describes  a  technique  that  has  been  developed  and  used  to 
rapidly  and  accurately  measure  these  properties  with  a  spatial  resolution  of 
1  mm.  Results  for  the  laser  synthesis  of  silicon  and  silicon  carbide 
powders  are  also  presented. 

III. 2.  EXPERIMENTAL  APPARATUS  AND  METHODOLOGY 

The  powder  synthesis  process  involves  passing  a  jet  of  reactant  gases 
through  a  CW  C02  laser  beam  inside  a  controlled  atmosphere  reaction  cell. 
Figure  III.l  is  a  schematic  of  the  reaction  zone.  The  reactant  gases  and 
the  C02  laser  beam  intersect  at  the  center  of  the  cell.  The  gas  is  heated 
rapidly,  and  reacts  to  form  small  particles,  which  are  swept  out  of  the  cell 
and  captured  in  a  filter.  Argon  is  injected  in  a  concentric  annular  flow  at 
the  base  of  the  nozzle  to  constrain  the  reactant  gases  in  a  narrow  stream 
and  to  carry  the  products  to  the  filter. 

By  this  method,  silane  gas  is  converted  into  silicon  powder  by  the 
simple  overall  reaction: 

SiH^g)  -  2H2(g)  +  Si(s)  . 


«5* 

«C- 


»t 


III.l.  Schematic  of  the  reaction  zone. 


The  resulting  silicon  powder  is  ideal  for  forming  reaction-bonded  silicon 
nitride  because  of  its  small  particle  size  (10-250  nm  diameter  depending  on 
process  parameters),  very  narrow  size  distribution,  high  purity,  and  loose 
agglomeration.  The  same  technique  has  been  used  to  produce  silicon  nitride! 
and  silicon  carbide  powders2. 

To  measure  the  process  parameters  of  interest,  an  apparatus  based  on  a 
HeN’e  probe  laser  was  developed.  The  extinction  and  scattering  of  the 
polarized  HeNe  beam  by  the  cloud  of  particles  provide  the  needed  data5,  by 
the  following  methodology. 

Optical  pyrometry  determines  the  brightness  temperature,  S,  of  a 
radiating  body,  at  a  single  wavelength.  The  absolute  emissivity,  e,  of  the 
area  being  viewed  must  be  known  before  the  true  temperature  can  be 
calculated.  This  can  be  readily  seen  from  Wien's  approximation  to  the  Plank 
Radiation  Law  written4  in  terms  of  the  radiance,  K.  ,  and  the  brightness 
temperature,  S,  and  the  true  temperature,  T: 


N, 
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If  the  emissivity  is  known,  the  true  temperature  T  can  be  found  by 
rearranging  Equation  1  to  yield: 


1  1  .  Xln  e 

-  =  —  +  - 

T  S  C2 


(2) 


The  accuracy  of  the  true  temperature  calculation  Is  essentially  equal  to  the 
accuracy  of  the  brightness  temperature  measurement  (t  10°C).  This  is 
because  the  emissivity  calculations  are  Insensitive  to  the  only  quantity 
that  is  not  directly  measured,  the  index  of  refraction  of  the  particles. 

The  emissivity  of  a  body  is  given  by  its  absorptivity, 

e  =a=l-x-p  (3) 

which  is  the  fraction  of  incident  light  that  is  neither  reflected,  p,  nor 
transmitted,  x.  For  a  diffuse  cloud  of  particles  such  as  exists  in  the 
reactions  of  interest,  light  is  scattered  in  many  directions  rather  than 
reflected  as  one  wave.  It  is  also  scattered  from  the  interior  of  the  cloud, 
rather  than  only  at  the  surface.  Thus,  the  emissivity  of  a  cloud  is:5 
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where  C  and  C  are,  respectively,  the  extinction  and  absorption  cross 

GXt  3  D  S 

sections  per  particle.  The  integrated  scattering  cross  section  C  is 

sea 

related  to  C  .  and  C  bv 
abs  ext  ' 

C  =  C  -  C  v  .  (5) 

sea  ext  abs 

The  cross  sections  are  functions  of  the  wavelength,  the  particle  size  and 
the  complex  index  of  refraction,  m  =  n  -  ik,  of  the  particles,  and  can  be 
calculated  using  the  Lorenz-Mie  scattering  equations.6  C  ^  re^ated  to 
the  transmissivity  through 

t  =  4:  =  exp(_ih'Cext)  (6) 

where  I0  and  I  are  the  incident  and  transmitted  intensities,  N  is  the 
particle  number  density,  and  1  is  the  thickness  of  the  cloud.  I'sing  these 
substitutions,  Equation  4  can  be  rewritten  as 

-  e  =  (1  11-4^)  * 

ext 

The  first  term  is  determined  from  a  transmission  measurement.  The  secom.-.' 
term  requires  that  the  scatter-extinction  ratio  be  measured. 

If  the  incident  and  scattered  light  are  plane  polarized  orthogonal  t 
the  scattering  plane,  the  scatter— extinction  ratio  is  defined  by: 

£  (?)  A2  i  (e) 

- i - = - i —  ,  (?■ 

Eofil  ln(Io/I)  4tt2  C 

ext 

where  the  quantities  on  the  left  side  are  measurable,  and  those  on  the  righ 

are  calculable.  £^(0)  is  the  power  scattered  at  angle  9  into  solid  angle 

for  an  incident  power  of  E0 ,  and  i^(@)  is  the  Mie  scattering  function  for 

light  scattered  at  angle  0.  Equation  8  is  written  for  the  case  where  the 

transmission  and  scattering  measurements  are  made  simultaneously,  so  that 

the  path  length  and  number  density  dependencies  cancel.  The  Mie  scattering 

function,  i^(0),  depends  on  the  particle  radius,  the  particle  index  of 

refraction,  and  the  wavelength  of  incident  light.6 

Values  for  C  ,  C  ,  and  i  (0)  are  calculated  as  a  function  of 
ext  sea  1 

particle  radius  for  the  values  of  0  and  A  used  in  the  experiment  with  an 
appropriate  index  of  refraction.  The  results  of  this  calculation  are 
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A. 


plotted  as  C  /C  versus  >.2  l,(0)/4r2  c  .  The  value  of  C  1C 

sea  ext  -l  ext  sea  ext 

corresponding  to  an  experimentally  measured  scatter-extinction  ratio  Is  then 

determined  from  this  plot.  The  emisslvlty  of  the  portion  of  the  cloud  on 

which  the  transmission  and  scattering  measurements  were  made  can  then  be 

calculated  from  Equation  7.  The  particle  si2e  can  be  determined  from  a  plot 

of  the  scatter-extinction  ratio  versus  particle  size.  The  cross  sections 

C  and  C  can  be  found  uslr®  similar  plots.  Once  C  and  the  width  of 
ext  sea  ext 

the  cloud  have  beer,  determine!,  the  particle  number  density  can  be 
calculated  using  Equation  . 

The  Lorenz-hle  scattering  theory  can  be  applied  to  the  scattering  of 
light  by  a  cloud  of  spherical  particles  If  the  following  three  conditions 
are  met:' 

(a!  Each,  particle  «ci?  'er<  separately;  that  is,  the  particles  are  not 
agglomerated,  and  are  serarat*1  at  least  3  radii;* 

(b)  There  is  no  optical  Interference  among  different  scattered  waves; 

(cl  There  Is  no  multiple  scattering. 

It  will  be  shown  In.  Section  III.'  that  condition  (a)  Is  true.  Condition  (b! 
Is  almost  always  valid  for  rando-  syste-s.’  But,  condition  (c)  requires 
that  the  turbidity  be  less  than  0.1,*  and  is  not  generally  met  by  regions  of 
the  clouds  of  particles  under  study.  The  most  significant  result  of  high 
turbidity  is  a  reduction  of  the  scattered  power  E  (?)  because  both  the  probe 
beam  and  the  scattered  light  from  the  interior  of  the  cloud  are  attenuated. 
For  small  angle  scattering,  this  effect  reduces  the  scattered  power  by  a 
factor  of  1  -  (I/I0)2/2  ln(I0/I).  (See  Appendix  III. A  for  the  derivation  of 
this  correction.)  Therefore,  using  2  fe^(e)  ln(I0/I)/[l  -  (X/X0)21  in  place 
of  £^(6)  in  Equation  8  allows  the  use  of  the  Mie  scattering  theory  even  in 
the  presence  of  high  turbidity. 

The  optical  properties  of  silicon  have  been  studied  extensively.  In 
the  visible  region  of  the  spectrum,  the  complex  index  of  refraction 
(m  ■  n  -  ik)  varies  strongly  with  wavelength,  temperature  and  degree  of 
crystallinity.8-10  Since  the  powders  being  studied  are  changing  temperature 
rapidly,  and  perhaps  crystallizing  as  they  grow,  a  series  of  calculations 
were  performed  using  various  possible  values  of  m. 

Figures  III, 2  and  III. 3  were  generated  from  computer  calculations1*1  of 
the  Mie  scattering  functions  and  cross  sections  using  three  different  values 
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III. 2. 


Scatter-extinction  ratio  versus 


Csca/Cexf 


III. 3,  Scatter-extinction  ratio  versus  particle  radius. 
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for  the  index  of  refraction.  An  index  of  n  =  4.1  -  0.1131  is  appropriate 

for  polycrystalline  silicon  at  500°C,1®  and  m  =  4.38  -0.99i  is  for  amorphous 

silicon  at  500°C.1°  An  extrapolated  value  of  n  =  4.38  -  0.26i  should  be 

valid  for  polycrystalline  silicon  around  1100°C  or  for  a  polycrystalline- 

amorphous  mix  at  some  lower  temperature.  Figure  III. 2  shows  the  calculated 

scatter-extinction  ratios  versus  C  /C  for  6  =  20°.  It  is  apparent  that 

sea  ext 

the  relationship  between  the  scatter-extinction  ratio  and  C  /C  is 

sea  ext 

independent  of  the  index  up  to  a  ratio  of  approximately  0.025.  The 

C  /C  values  remain  close  to  one  another  for  the  two  polycrystalline 
sea  ext 

indices  for  higher  values  of  the  ratio.  At  a  scatter-extinction  ratio  of 
0.04,  the  maximum  uncertainty  in  e  results  in  a  true  temperature  difference 
of  only  26°C  for  a  brightness  temperature  of  1100°C  (~  2.5TO.  This 
indicates  that  even  though  the  index  as  a  function  of  position  in  the 
reaction  zone  is  not  precisely  known,  the  emissivity  and  the  temperature  can 
still  be  determined  to  high  accuracy. 

Figure  III. 3  is  a  plot  of  the  scatter-extinction  ratio  versus  particle 
radius  for  the  same  three  indices.  A  much  larger  dependence  on  index  is 
noted.  Although  this  introduces  some  uncertainty  to  the  absolute  particle 
size  and  number  density  calculations,  this  is  not  a  serious  problem  since 
the  final  particle  size  can  be  measured  by  BET  surface  area  analysis  and 
transmission  electron  microscopy.  Without  question,  the  scatter-extinctior. 
measurements  can  be  used  to  monitor  changes  in  particle  size  and  number 
density  throughout  the  reaction.  It  may  also  be  possible  to  extract 
information  on  the  variation  of  the  Index  of  refraction  of  the  particles 
within  the  reaction  zone  from  the  scatter-extinction  measurements  as  a 
function  of  position. 

III. 3.  EXPERIMENTAL  PROCEDURE 

A  HeNe  laser  is  used  as  the  probe  in  an  apparatus  that  simultaneously 
measures  the  extinction  of  the  HeNe  laser  beam  and  the  intensity  of  the 
scattered  light  at  an  angle  of  20°.  Figure  III. 4  Is  a  schematic  of  the 
optical  layout. 

The  5  mW  vertically  polarized  HeNe  laser  is  mounted  on  a  cathetoneter, 
enabling  measurements  to  be  made  as  a  function  of  height  by  simply 
translating  the  laser.  The  transmitted  HeNe  laser  light,  I,  Is  measured  by 

I 
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Scatter-extinction  optical  layout. 


a  silicon  photodiode  nounted  on  an  x-y  translator.  The  scattered  HeN'e 
light,  E^(20°),  is  detected  by  an  RCA  C53050  photomultiplier.  An  8.13-mn- 
diameter  aperture  determines  the  solid  angle  viewed,  which  is  2.23  x  10~4 
steradians.  Focussing  optics  are  not  used  to  minimize  alignment  problems  as 
the  probe  laser  is  moved.  A  632.8  nm  bandpass  filter  passes  the  HeKe  laser 
light,  but  rejects  most  of  the  thermal  emission  from  the  reaction  2one.  A 
vertical  polarizer  analyzes  the  polarization  of  the  scattered  light.  The 
incident  probe  laser  power,  E0,  was  measured  by  directing  a  calibrated 
fraction  of  the  HeNe  laser  beam  through  the  same  collection  filters  to  the 
photomultiplier  using  a  bean  splitter  temporarily  mounted  at  the  center  of 
the  reaction  cell. 

A  micro  optical  pyrometer  (made  by  Pyrometer  Instruments,  Inc.)  Is 
mounted  alongside  the  photodiode  so  that  it  views  the  reaction  zone  in  the 
scattering  plane.  The  field  of  view  of  the  pyrometer  is  such  that  a  1  mm2 
area  is  monitored.' 

At  each  height  the  following  sequence  of  adjustments'  and  neasurener. t? 
is  required:  the  cathetometer  is  set  to  the  desired  height;  the  photodiode 
is  adjusted  with  its  x-y  translator  for  peak  signal;  I  and  f^(2n°)  art- 
measured  with  the  lock-in  amplifier;  the  optical  pyrometer  is  focussed  on 
the  point  where  the  HeN'e  beam  intersects  the  reaction  zone,  then  the  HeN'e 
laser  is  blocked  while  a  temperature  reading  is  made.  In  this  way  the 
scattering,  extinction,  and  brightness  temperature  measurements  are  all  made 
on  the  identical  volume  of  the  reaction  zone. 

III. 4.  RESULTS  AND  DISCUSSION 

Measurements  were  made  on  laser-heated  silicon  powder  synthesis 
reactions  under  a  range  of  process  conditions.  Table  III.l  lists  the 
parameters  of  the  reactions  studied  and  the  equivalent  particle  sizes 
determined  by  BET  surface  area  analysis  for  similar  runs.1 

Figure  III. 5  is  the  measured  scattered  power,  £^(20°),  and  transmitted 
intensity,  I,  in  millivolts  RMS  for  Run  631S  as  a  function  of  position  above 
the  reactant  gas  nozzle.  The  scatter-extinction  ratio  is  also  plotted. 
Figure  III. 6  shows  the  emissivity,  t,  the  brightness  temperature,  S,  and  the 
true  temperature,  T,  as  a  function  of  position  for  the  same  reaction. 


Temperature  (°C)  n'i,liv°"*  rms 


£^(20°),  and  scatter-extinction  ratio  for  run  631S 


Distance  from  nozzle  (mm) 


III. 6.  Temperature  and  emissivity  for  run  631S. 


Temperature  (°C) 


Temperature  and 


emissivi ty  for  run  654S. 


Temperature  versus  distance  from  nozzle. 
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III. 9 


.  Particle  size  versus  distance  from  nozzle. 


III. 10 


Particle  number  density  versus  distance  from  nozzle 


Figure  III. 7  shows  c,  S,  and  T  for  Run  654S.  These  plots  are  representative 
of  the  data  from  the  other  runs.  Figure  III. 8  summarizes  the  true 
temperature  distributions  for  all  of  the  reactions  studied.  Figure  III. 9 
summarizes  the  calculated  particle  sizes. 

The  number  density  of  particles  as  a  function  of  position  in  the 
reaction  zone  is  presented  in  Figure  III. 10  for  Run  654S.  Two  different 
calculations  are  presented,  since  the  calculated  number  density  is  a  strong 
function  of  the  value  of  the  index  of  refraction  used. 

Figures  III. 8  and  III. 9  reveal  that  the  heating  of  the  reactants  is 
rapid,  and  the  reaction  is  often  underway  before  the  reactants  intersect  the 
C02  laser  beam.  In  the  case  of  Runs  650S  and  631S,  the  measured 
temperatures  at  the  edge  of  the  laser  beam  are  over  900°C. 

Possible  heating  mechanisms  of  the  gas  below  the  CC>2  laser  beam  are 
thermal  conduction  from  the  hot  zone,  and  resonant  emission  from 
vibrationally  excited  Si H4  molecules  in  the  hot  zone.  The  maximum  rate  of 
infrared  emission  from  the  reaction  zone  can  be  calculated  from  the  Planck 
radiation  law.  At  X  =  4  ym,  and  1100°C,  the  emission  rate  is  about 
3  W/cm2/ym.  If  the  width  of  the  4  ym  silane  absorption  band  is  0.33  yr, 

1.0  W/cr.2  could  be  transferred  by  radiation.  Hot  bands,  and  the  10  ym  band 
will  also  contribute  to  some  extent.  Thermal  conductivity  across  a  3-mr  gap 
of  silane  is  about  2  W/cm2  for  a  1000°  temperature  difference.  From  this 
simple  analysis  it  appears  that  both  thermal  conductivity  and  vibrational 
radiation  from  excited  silane  molecules  are  significant  In  the  preheating  of 
the  reactants. 

Figure  III. 8  also  indicates  that  for  all  reaction  conditions  studied 
the  temperature  of  the  silicon  powder  begins  decreasing  while  the  particles 
are  still  In  the  C02  laser  beam,  suggesting  that  the  particles  themselves 
are  not  strongly  heated  by  the  laser.  As  Figure  III. 9  shows,  virtually  all 
of  the  particle  growth  takes  place  within  the  C02  laser  beam.  Note  that  the 
continuing  growth  in  Run  630S  Is  due  to  the  high  silane  concentration  and 
flow  rate.  Once  most  or  all  of  the  silane  has  reacted  to  form  solid 
silicon,  the  gas/powder  mixture  stops  absorbing  and  starts  to  cool  off. 

Since  silane  decomposition  starts  around  800°C,  most  of  It  is  consumed  at 
temperatures  significantly  below  the  silicon  melting  point. 

Since  the  particle  size  determinations  are  sensitive  to  the  index  of 
refraction,  the  quantitative  conclusions  based  on  these  measurements  must  be 
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more  tentative.  It  is  clear  that  the  particles  grow  rapidly,  and  that  the 
growth  process  usually  terminates  within  the  CC>2  laser  beam  volume. 
Significant  growth  past  the  laser  beam  only  occurs  at  the  high  concentration 
and/or  flow  rates  (as  in  Run  630S  and  631S).  F<-  r  the  undiluted  0.2  atm 

runs,  agreement  between  these  measurements  and  BET  size  measurements  is 
good,  suggesting  that  the  polycrystalline  indices  extrapolated  above  SOO^C 
are  the  best  choice.  For  diluted  reactions,  which  are  not  as  hot,  better 
agreement  is  obtained  when  the  lower  temperature  indices  are  used.  The 
overall  good  agreement  indicates  that  the  high  turbidity  correction  is 
reasonably  good  for  turbidities  up  to  about  2. 

The  run  at  0.6  atm  (650S)  requires  additional  interpretation.  The 
measured  particle  size  is  much  smaller  than  the  BET  results.  However,  since 
the  turbidity  of  this  reaction  is  around  6,  it  is  quite  possible  that  the 
turbidity  correction  is  not  accurate  at  very  high  turbidities.  The  other 
two  conditions  for  the  validity  of  Mie  scattering  should  still  hold. 

Another  explanation  is  that,  for  such  a  turbid  cloud,  virtually  all  of  the 
detected  scattering  comes  from  the  outer  edge,  where  the  silane  reactant  is 
the  most  strongly  diluted  with  argon.  Since  diluted  conditions  produced 
smaller  particles,  it  is  likely  that  small  particles  are  created  at  the 
edges  of  the  reaction  zone.  The  effect  of  this  would  be  to  reduce  the 
scatter-extinction  ratio,  leading  to  a  small  calculated  particle  size. 

Finally,  the  leveling  off  of  the  scatter  extinction  ratio,  hence  the 
particle  size,  under  all  conditions  indicates  that  the  large  loose 
agglomerates  that  are  typical  of  the  powder  after  it  has  been  collected  in 
the  filter  are  produced  by  the  collection  process.  This  justifies  the 
assumption  that  condition  (a)  for  Mie  scattering  is  valid.  It  also  suggests 
that  a  liquid  collection  process  should  produce  a  monodisperse  suspension  of 
silicon  particles.  This  would  allow  the  fabrication  of  uniform  high  density 
compacts,  which  would  have  much  improved  sintering  properties. 

The  particle  number  density  calculations  are  extremely  sensitive  to  the 
value  of  the  refractive  index  used  in  the  calculation.  A  calculation  based 
on  a  constant  index  Indicates  that  the  number  density  decreases  from  a 
maximum  value  at  the  inception  of  the  reaction  to  a  constant  value 
(Figure  III. 10).  Several  factors  can  contribute  to  the  apparent  decrease  in 
number  density  from  levels  calculated  at  the  base  of  the  reaction  zone. 
Although  expansion  of  the  reactant  gases  with  heating  is  a  factor,  the 


number  density  change  is  too  great  for  this  to  be  solely  responsible.  Small 
particles  could  be  agglomerating  as  they  move  into  the  C02  laser  bear.  We 
feel  this  is  unlikely  because  there  is  no  reason  why  the  agglomeration 
process  would  terminate  abruptly,  particularly  before  the  reaction  had 
neared  completion.  Figure  III. 9  shows  that  the  major  portion  of  the 
reaction  occurs  via  growth  after  the  number  density  has  reached  a  constant 
value.  The  final  explanation  is  that  the  particles  are  amorphous  at  the 
base  of  the  reaction  zone  and  crystallize  as  they  grow  while  travelling  into 
the  hotter  regions  of  the  zone.  Under  this  assumption,  the  number  density 
curve  would  be  a  combination  of  the  amorphous  curve  and  the  polycrystalline 
curve,  as  shown  by  the  dashed  line  in  Figure  III. 10.  The  transition 
corresponds  to  the  region  within  the  zone  where  the  amorphous  particles 
crystallize.  The  temperature  at  the  end  of  the  transition,  4.5  mm  from  the 
nozzle,  was  1120°C.  At  the  beginning  point,  3  mm  from  the  nozzle,  the 
reaction  zone  was  not  luminous,  so  the  temperature  there  is  unknown.  A 
linear  extrapolation  of  the  temperature  indicates  750°C,  which  is  probably 
an  upper  limit.  This  low  temperature  is  consistent  with  the  formation  of 
amorphous  silicon,  just  as  1120°C  is  consistent  with  polvcrystalline  silicon 
formation.  The  formation  of  amorphous  powders  and  their  progressive 
crystallization  is  therefore  probably  responsible  for  the  apparent  number 
density  decrease  with  height  when  a  constant  refractive  index  is  assured. 

The  dashed  line  in  Figure  HI. 10  indicates  a  relatively  constant  number 
density. 

Preliminary  scatter-extinction  and  temperature  measurements  have  been 
performed  on  SiC  powders  made  from  SiH^  and  CjH^.  Figu~e  III. 11  shows  the 
brightness  temperature,  true  temperature,  and  emissivity  for  a  typical 
silicon  carbide  synthesis  run.  The  run  conditions  are  listed  in 
Table  III. 1 .  As  was  observed  In  the  silane  reaction,  the  heating  of  the 
reactants  is  very  rapid  and  the  reaction  is  initiated  significantly  before 
the  reactants  enter  the  laser  beam.  The  emissivity  increased  rapidly  until 
a  nearly  constant  value  was  achieved.  The  temperature  distribution  within 
the  laser-irradiated  region  was  sufficiently  similar  in  shape  to  the 
Gaussian  intensity  distribution  of  the  laser  beam  to  suggest  that,  unlike 
the  silane  synthesis  process,  the  SiC  particles  are  heated  directly  by  the 


III. 11.  Temperature  and  emissivity  for  Run  #  910S1C. 


III. 12.  Particle  size  and  number  density  for  Run  #  910S1C 


The  calculated  particle  radius  is  shown  in  Figure  III. 12.  The  hi<-h 
temperature  optical  properties  of  SiC^2  are  not  precisely  known,  and  the 
imaginary  part  varies  strongly  with  small  departures  from  stoichiometry. 
Therefore  a  series  of  calculations  using  various  values  for  the  refractive 
index  was  carried  out;  the  set  of  calculations  that  produces  a  final 
particle  size  that  agrees  the  closest  with  the  BET  value  is  plotted  in 
Figure  III. 12.  It  is  evident  that  the  nucleation  and  growth  processes  are 
nearly  complete  by  the  time  the  particles  reach  the  laser  beam,  suggesting  a 
very  rapid  growth  rate.  Better  data  on  the  SiC  optical  properties  should 
enable  more  accurate  comparisons  between  the  scatter-extinction  results  and 
the  particle  size  determined  by  other  analytical  techniques. 

Figure  III. 12  also  shows  a  plot  of  the  particle  number  density,  N ,  as  a 
function  of  position  In  the  reaction  zone.  The  calculated  number  density 
decreased  from  a  maximum  value  at  the  base  of  the  reaction  zone  and 
continues  to-  decrease  slightly  as  the  particles  travel  further  fror  the 
inlet  nozzle.  The  decrease  in  particle  number  density  may  indicate  the 
particles  are  agglomerating  within  the  reaction  zone  or,  as  suggested  in  thc- 
silicon  synthesis,  the  particles  may  initially  form  as  amorphous  SiC  and 
crystallize  in  the  hotter  regions  of  the  reaction  zone.  A  transition  from 
silicon  rich  to  nearly  stoichiometric  SiC  would  also  explain  this  decrease. 

A  more  detailed  study  is  required  to  interpret  these  results  with  respect  to 
the  particle  formation  mechanisms. 

III. 5.  CONCLUSION 

The  analytical  techniques  presented  here  are  applicable  to  the  accurate 
determination  of  the  temperature  distribution  and  average  particle  size  and 
number  density  of  a  cloud  of  particulates,  including  silicon,  silicon 
carbide  or  silicon  nitride,  produced  in  a  C02  laser-driven  reaction.  While 
the  technique  has  been  primarily  used  to  study  the  production  of  Si  powder 
from  silane  gas,  an  equivalent  study  of  the  SiC  and  SijN^  reactions  requires 
only  some  idea  of  their  high-temperature  optical  properties.  This  technique 
is  relatively  insensitive  to  the  optical  thickness  of  the  cloud,  or  to 
whether  the  individual  particles  primarily  absorb  or  scatter  light  (although 
a  finite  absorption  cross  section  is  required). 
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APPENDIX  III. A. 


Correction  for  High  Turbidity  In  Scattering  Measurements 

The  derivations  from  Lorenz-Mie  Scattering  theory  are  valid  only  when 
the  three  conditions  cited  in  the  text  for  single  scattering  are  met.  If 
the  probe  beam  is  attenuated  significantly  by  the  cloud,  the  last  volume 
element  will  experience  an  incident  power  less  than  E0,  so  it  will  scatter 
less  light  than  does  the  first  volume  element.  In  small  angle  scattering, 
the  light  from  the  last  volume  element  of  the  cloud  will  have  to  travel  back 
through  the  cloud,  being  further  attenuated.  From  Lorenz-Mie  theory,  the 
light  scattered  at  position  p,  t(G,p)  is 

t(9,p)  =  N’Csca  ( 6 )  E(p)  fiC  ip  ,  (A-i) 

where  N  is  the  particle  number  densitv,  C  (©)  is  the  scattering  cross 
f  '  sea 

section  at  angle  © ,  E(p)  is  the  probe  power  at  position  p,  and  fifi  and  6p  are 
the  solid-angle  and  volume  depth  respectively.  If  the  cloud  is  optically 
thin,  E(p)  =  Ec,  and  the  scattered  light  is  not  attenuated  as  it  leaves  the 
cloud,  so  the  total  scattered  light  from  the  entire  cloud  of  thickness  E  is 

D 

£(e)  =  f  E(e,p)  dp  -  n  c  (e)  e0  fin  d  .  (a-2) 

0 '  sca 

where  D  is  the  thickness  of  the  cloud.  If  the  cloud  is  not  optically  thin, 
then  the  power  at  position  p  is 

-NC  p 

E(p)  =  E0  e  eXt  (A-3) 

-NC  L 

6X  t 

and  the  scattered  light  is  attenuated  by  e  ,  where  L  is  the  distance 

out  of  the  cloud  at  an  angle  0  from  position  p.  Figure  III.A-1  illustrates 

these  relations  for  a  circular  cloud  with  the  probe  beam  passing  along  the 

diameter,  D.  The  total  scattered  power  is 

D  D  -NC  p  -NC  L 

£(e)  =  /  £<e,p)dP  =  |nc  (e)  fin  E„e  ext  e  ext  dp  .  (a-4) 

0  0  sca 

For  6  *  20°,  L  is  equal  to  p  within  5%  so  Equation  A-4  can  be  integrated: 


46 


i 


E(6)  =  N  C  (e)  fin  Ec  '  e 
sea  0 


-2NC 


1-e 


ex 


N  c  (e)  fin  EeD  f— - =- 

sca  v  2NC  D 

ext 


From  the  Beer-Lambert  relation: 

I 

I« 


-NC  D 

=  e  ext 


or 


ln(Io/I)  =  NC  D 
ext 


Equation  A- 5  can  therefore  be  rewritten  as 

1  -  (I/l0)2 

£(e)  =  n  c  (e)  fin  e0d  (— • - — t- — 

sea  1  2  In  (I0/I) 


III.A-1.  Scattering  diagram  with  cylindrical  symmetry 


For  the  case  where  the  incident  and  scattered  waves  are  orthogonal  to  the 

scattering  plane,  — —  i,  (6)  is  substituted  for  C  (0),  leading  to 

41,2  1  Sca 


2^(8)  ln(Io/I )  X2  ii(0)  N  D 

-  —  -  ■-  m  . -  ■  —  ■  — 

Ee6n(l  -  (I/Io)2)  4 1,2 


(A-9) 


Again,  the  right  side  is  calculable  from  the  Mie  equations,  the  left  is 

measured.  Dividing  by  ln(I0/I)  *  NC  D  leads  to  the  scatter-extinction 

ext 

ratio  corrected  for  high  turbidity, 


2  £^(e)(measured) 


X2  ii(e) 

4 1,2  C  . 
ext 


E„6fl  (l-(I/Io)2) 


(A-10) 


CHAPTER  XV 


HOMOGENEOUS  HU CLEAT ION  AND  GROWTH  OP  SILICON  POWDER 

PROM  LASER  HEATED  SiH, 

A 

by 

R.  A.  Marra,  J.  H.  Plint  and  J.  S.  Haggerty 

ABSTRACT 

Nucleation  and  growth  kinetics  of  Si  particles  from  C02  laser  heated 
S1H4  gas  have  been  measured  by  means  of  a  HeNe  laser  light  scattering  and 
transmission  technique.  This  has  permitted  the  reaction  kinetics  to  be 
monitored  in  terms  of  temperature  and  reactant  partial  pressures. 

With  heating  at  rates  between  105-107°C/s  silicon  particles  nucleate 
homogeneously  at  rates  in  excess  of  lO^/cm^s.  The  rates  are  in  reasonable 
agreement  with  Volmer-Weber  theory.  The  rate  is  controlled  by  the 
concentration  of  silicon  atoms  in  the  vapor  phase;  the  high  level  of 
supersaturation  results  in  a  negligible  barrier.  Nucleation  ceased  after  the 
reaction  proceeded  only  a  few  percent  apparently  resulting  from  the  rapid 
decrease  in  supersaturation  and  relatively  favored  heterogeneous  growth. 
Growth  rate  is  controlled  by  vapor  phase  transport  of  reactants  to  the 
particle  surfaces.  The  condensation  coefficient  is  estimated  to  be  0.002- 
0.02.  The  growth  process  is  terminated  by  the  depletion  of  reactants;  the 
final  particle  size  is  established  by  the  initial  SiH^  concentration  and  the 
number  of  nuclei  formed. 

Separation  of  the  nucleation  and  growth  processes  permits  final 
particle  dimensions  to  be  controlled  through  process  variables.  It  also 
results  in  uniform  particle  dimensions  needed  for  ceramic  processes  based  on 
ordered  dispersions. 

IT. 1.  INTRODUCTION 

Reconsideration  of  the  relationships  between  the  properties  of  fired 
ceramic  bodies  and  the  characteristics  of  their  constituent  powders  suggested 
that  uniform  diameter  powders  should  be  substituted  for  the  traditionally 
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used  powders  with  broad  particle  size  distributions.1  Achieving  reliable 
slip  stability  also  requires  highly  controlled  powder  purity,  stoichiometry, 
crystallinity,  agglomeration,  shape  and  size.  Because  comminution  will  not 
produce  submicron  powders  exhibiting  all  of  these  characteristics,  they  must 
be  synthesized  In  a  manner  that  achieves  these  properties  simultaneously.  In 
turn,  this  requires  that  nucleatlon,  growth,  and  agglomeration  processes  be 
understood  and  controlled. 

2,3 

Gas  phase  synthesis  generally  achieves  the  required  attributes.  The 

establishment  of  high  supersaturation  levels  results  In  small,  equiaxial 
particles  because  nucleatlon  rates  are  high  and  growth  is  Isotropic.4 
Agglomeration  can  be  controlled  with  dilution  and  by  minimizing  exposure  to 
high  temperatures.  Overall  purities  and  compositions  are  controlled  through 
the  reactants. 

5  6 

Numerous  techniques  are  utilized  to  synthesize  powders  from  vapor 
phase  reactants.  These  include  tube  furnaces,  arc  and  rf  plasmas,  flash 
thermolysis,  evaporation/condensation  in  inert  atmospheres  and  shock  tubes. 
While  all  achieve  some  of  the  required  process  attributes,  each  but  the  last 
falls  to  achieve  them  all.  Usually,  the  thermal  profiles  within  reaction 
zones  result  in  nonuniform  nucleatlon  and  growth  histories,  or  particles  do 
not  cool  rapidly  enough  to  avoid  formation  of  hard  agglomerates. 

The  laser-driven  gas  phase  reaction  process2  permits  the  achievement  of 
all  the  desired  powder  characteristics. 3  In  this  process,  an  optical  energy 
source,  a  C02  laser,  Is  used  to  Initiate  and  sustain  a  vapor  phase  chemical 
reaction.  The  gaseous  molecules  absorb  the  laser  energy,  which  Is 
thermallzed  through  molecular  collisions.  The  reaction  vessel  walls  remain 
cold,  avoiding  contamination  and  heterogeneous  nucleatlon  sites.  The 
reaction  zone  Is  well  defined  by  the  Intersection  of  the  laser  beam  and 
reactant  gas  stream.  A  wide  range  of  uniform  and  precise  heating  rates  with 
controlled  maximum  temperatures  can  be  achieved  by  varying  the  laser 
intensity  as  well  as  the  gas  mixture,  pressure,  and  flow  rates.  The  amount 
of  energy  supplied  to  the  gas  and  the  resulting  heating  are  defined  by  the 
absorptivity  of  the  gaseous  products.  If  the  wavelength  of  the  emitted  light 
and  the  absorption  lines  of  the  gaseous  molecules  are  In  near  coincidence, 
the  optical- to- thermal  efficiency  Is  very  high.  Importantly,  the  physical 
geometry  of  the  laser  synthesis  process  enables  nucleatlon  and  growth 


kinetics  to  be  monitored  in  terms  of  measured  temperatures  and  calculated 
reactant  partial  pressures. 

The  research  described  in  this  paper  involves  the  nucleation  and  growth 

mechanisms  of  silicon  powders  synthesized  from  laser-heated  silane  gas.  The 

pyrolysis  of  silane  to  form  silicon  was  modeled  first  because  the  reaction 

path  is  relatively  simple  compared  to  SijN^  and  SiC  synthesis  processes,  and 

because  of  the  technical  importance  of  high-purity  silicon  powder  for 

2  7 

reaction- bonded  silicon  nitride.  The  empirical  modeling  *  of  the  synthesis 

process  and  the  effect  process  variables  have  on  the  reaction  zone  and 

3  8  9 

particle  characteristics  ’  *  were  reported  in  earlier  publications.  The 
experimental  techniques  used  to  measure  particle  diameter,  number  density  and 
temperature  were  developed  by  Flint^0  and  are  reported  in  detail 
elsewhere. 1 1 

IV. 2.  EXPERIMENTAL  AND  ANALYTICAL  PROCEDURES 

The  powder  synthesis  cell  is  shown  in  Figure  IV. 1.  A  cross-flow 
configuration  is  used  in  which  a  laser  beam  orthogonally  intersects  the 
reactant  gas  stream.  The  gas  stream  is  made  optically  thin  so  all  portions 
of  the  gas  stream  are  uniformly  exposed.  The  laser  beam  enters  and  leaves 
the  cell  through  KC1  windows.  The  silane  reactant  gas,  under  some  conditions 
diluted  with  argon,  enters  the  cell  through  a  1.5  mm  ID  stainless  steel 
nozzle  located  6  mm  from  the  beam  center  line.  The  powders  produced  in  the 
reaction  zone  are  carried  to  a  microfiber  filter  by  the  product  gases  and  a 
surrounding  coaxial  stream  of  argon. 

Synthesis  process  variables  were  manipulated  to  determine  the 
interactions  between  the  reaction  zone  and  particle  characteristics.  The 
synthesis  conditions  are  listed  in  Table  IV. 1. 

The  reaction  temperature  in  the  laser  powder  process  is  not  an 
extrin8lcally  controllable  parameter  but  is  dependent  on  the  process 
variables.  Because  the  temperature  within  the  reaction  zone  determines  the 
thermodynamic  driving  forces,  vapor  phase  reaction  kinetics,  and  particle 
crystallinity,  it  is  important  to  measure  the  temperature  distribution  within 
the  reaction  zone. 

A  micropyrometer  was  used  to  determine  the  brightness  temperature  of  the 
hot  particles  in  a  small  region  of  the  reaction  zone.  True  temperatures  were 
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calculated  using  the  local  emlsslvlty  of  the  particulate  cloud  determined  by 

probing  the  reaction  zone  with  a  HeNe  laser  and  measuring  the  scattered  and 

10.11 

transmitted  light  intensities.  The  ratio  of  the  scattering  to 

extinction  cross-sections  determined  from  experimental  results  and  Lorenz-Mle 
scattering  theory  permits  the  emlsslvlty  of  the  particulate  cloud  to  be 
calcula ted. 

The  scatter-extinction  method  also  permits  the  particle  number  density 

and  particle  diameter  to  be  measured  as  functions  of  position  In  the  reaction 
10  11 

zone.  The  Mie  scattering  function  and  extinction  cross-sections  are 

highly  dependent  on  the  particles'  refractive  index.12  Calculated  and 

measured  values  of  final  particle  diameters  were  compared  to  insure  that 

8.1011 

correct  refractive  indices  were  used  in  calculations.  Local  reactant 

partial  pressures  were  determined  by  the  difference  between  the  Initial 
reactant  concentration  and  the  local  product  concentration  calculated  from 
the  local  particle  number  density  and  diameter.  Combined  with  calculated  gas 
velocities,  this  Information  permits  nucleation  and  growth  kinetics  to  be 
determined  as  functions  of  temperature  and  reactant  partial  pressures. 

IV.  3.  RESULTS  AND  OBSERVATIONS 

IV. 3. A.  Temperature  Measurements 

The  general  features  of  the  temperature  distribution  In  the  reaction 

zones  were  similar  for  all  synthesis  conditions.  Figure  IV. 2  shows  the 

emlsslvlty,  brightness  temperature,  and  true  temperature  for  Run  631S  (see 

Table  IV. 1  for  synthesis  conditions).  For  this  and  some  other  run 

conditions.  It  is  evident  that  the  reaction  was  well  underway  before  the 

8  10 

reactants  reached  the  CO^  laser  beam.  Thermal  analyses  *  showed  that  at 
low  flow  rates,  the  silane  gas  is  preheated  by  conduction  of  heat  toward  the 
inlet  nozzle.  The  temperature  gradient  Increases  rapidly  and  reaches  a 
maximum  before  the  gas  reaches  the  beam  center.  There  Is  no  apparent 
correlation  between  the  particle  temperature  and  the  laser  Intensity 
distribution,  Indicating  that  the  silicon  particles  are  not  strongly  heated 
by  the  C02  laser  beam;  this  is  expected  since  the  absorption  efficiency  of 
50  nm  silicon  particles  to  10.6  pm  light  Is  quite  low. 

The  time-temperature  history  conforms  to  our  requirements.  Heating 
rates  are  high,  typically  106-108°C/s.  The  particles  experience  no  sustained 
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Figure  IV.  2.-  True  temperature,  brightness  temperature,  and  emisslvity  as  a 
function  of  distance  above  the  inlet  nozzle  for  Run  631S. 


Figure  IV. 3.  Particle  size  calculated  from  the  scatter-extinction  results  as 
a  function  of  distance  above  the  inlet  nozzle  for  various 
synthesis  runs. 
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exposure  to  elevated  temperatures.  The  level  of  the  maximum  temperature  is 
low  compared  to  many  techniques  and  its  value  can  be  controlled  over  a  wide 
range  to  an  accuracy  of  t  20°C.  Finally,  the  cooling  rates  are  rapid  enough, 
~  105°C/s,  to  retain  metastable  structures. 

IV. 3. B.  In  Situ  Nucleation  and  Growth  Heaaurenents 

Figure  IV. 3  summarizes  the  particle  size  measurements  as  a  function  of 
distance  above  the  inlet  gas  nozzle.  The  particles  grow  at  a  rapid,  nearly 
constant  rate,  ~  106A/s,  to  their  final  diameter  while  still  in  the  laser 
beam.  Under  certain  conditions,  the  reaction  begins  prior  to  reaching  the 
laser  beam.  Comparison  of  the  results  in  Figure  IV. 3  with  corresponding 
temperature  profile  measurements  shows  that  the  maximum  temperature  is 
reached  at  a  point  corresponding  to  the  termination  of  the  growth  process. 

For  runs  using  undiluted  silane  reactant  gas  at  2  x  104Pa  (0.2  atm),  the 
closest  agreement10  between  the  directly  measured  TEM  and  BET  diameters  and 
the  calculated  diameters  resulted  when  the  particles  were  assumed  to  have  a 
refractive  index  of  m  =  4.38-0.261  suitable  for  polycrystalline  silicon15  at 
around  1100°C  (Runs  630S  and  631S).  For  runs  diluted  with  argon,  (Runs  654S 
and  Runs  634S),  the  best  agreement10  was  obtained  using  a  refractive  index 
corresponding  to  a  somewhat  lower  temperature;  both  were  consistent  with  the 
temperature  measurements  and  the  crystal  structures  of  the  final  particles.5 
For  the  run  performed  at  0.6  atm  (Run  650S),  the  calculated  particle  size  was 
smaller  than  the  measured  diameter  for  all  reasonable  refractive  index 
values.  For  these  process  conditions,  90%  of  the  incident  HeNe  laser  light 
was  attenuated  in  the  first  millimeter  of  the  particulate  plume;  thus,  the 
majority  of  the  scattered  probe  signal  originates  from  the  outer  edge  of  the 
particle  cloud.  This  outermost  region  is  most  highly  diluted  by  the 
surrounding  argon  gas  flow  and  is  subjected  to  a  lower  C02  laser  intensity 
than  is  the  reactant  gas  that  passes  through  the  beam  center.  Therefore,  it 
is  likely  that  the  gas  stream  temperature  in  this  outer  zone  and  resulting 
particle  diameter  will  be  significantly  smaller  than  the  average  values. 

Particle  number  density  calculations  employ  the  HeNe  path  length  and  the 
results  of  the  scatter-extinction  measurements.  Because  it  is  difficult  to 
measure  the  path  length  l  through  the  reaction  zone  accurately,  the  product 
of  the  number  density  and  the  path  length  was  calculated  and  is  reported. 

For  most  synthesis  conditions,  the  path  length  was  relatively  constant  and 
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approximately  equal  to  1  cm,  so  the  number  density  is  well  represented  by 
this  product. 

Figure  IV. 4  shows  the  calculated  number-density»path-length  product,  Nt, 
as  a  function  of  position  above  the  gas  inlet  nozzle  for  the  same  conditions 
as  in  Figure  IV. 3. 10 

The  calculated  Ni  values  decrease  from  maximum  values  at  the  initiation 
of  the  reaction  zone  until  they  reach  constant  values  further  into  the 
particulate  cloud.  The  volumetric  expansion  of  the  heated  gas  has  only  a 
very  small  effect  on  the  Ni  product.  The  decrease  in  Ni  might  indicate  that 
particles  coagulate  as  they  penetrate  the  reaction  zone;  however,  an  abrupt 
termination  of  the  agglomeration  process  is  unlikely.  Alternatively,  the 
number  density  remains  essentially  constant  if  the  particles  are  assumed  to 
have  a  refractive  index  appropriate  for  amorphous  silicon  when  they  form  and 
then  change  to  a  polycrystalline  index  with  further  penetration  into  the 
reaction  zone.  This  behavior  indicates  that  particles  form  as  amorphous 
silicon  and  then  crystallize  in  the  deeper  regions  of  the  reaction  zone;  this 
is  consistent  with  the  final  particle  microstructure  and  possible  particle 
formation  mechanisms.14 

The  constant  number  density  indicates  that  for  these  synthesis 
conditions  the  particles  do  not  agglomerate  within  the  observed  region  of 
reaction  zone.  These  results  also  indicate  that  nucleatlon  occured  only  at 
the  very  base  of  the  reaction  zone,  ceasing  before  significant  particle 
growth  occurred.  Both  results  are  critical  for  producing  powders  with  the 
desired  characteristics. 

IT. 4.  PARTICLE  RUCLEATION  AND  GROWTH  MECHANISMS 

IT. 4. A.  Hue lea t ion 

The  exact  reaction  path  and  rate-controlling  step  for  silane 
decomposition  are  not  known  with  certainty;  both  simple  decomposition1 5  and 

polymerization  mechanisms1 6  have  been  proposed.  Several  lnvestiga- 

17  18  19 

tions  of  C02  laser- induced  silane  pyrolysis  have  shown  the  absence  of 

higher  silanes  in  the  reaction  products,  suggesting  that  the  laser  Induced 
pyrolysis  occurs  by  the  simple  decomposition  reaction; 
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SiH4(g)  -  Si( s)  +  2H2(g)  .  (1) 

We  have  developed  our  analyses  on  this  basis.  Because  this  reaction  has  a 
negative  free  energy  and  is  exothermic, 20  it  is  apparent  that  kinetic  rather 
than  thermodynamic  factors  would  dominate  nucleation  and  growth  processes. 

For  nucleation  to  occur,  the  vapor  phase  must  be  sufficiently 
supersaturated  to  overcome  the  free  energy  barrier  created  by  the  new 
surface  of  a  nucleus.  The  nucleation  rate  depends  on  the  number  of  silicon 
atoms  in  the  vapor  phase,  since  nucleation  occurs  by  an  accretion  process 
limited  by  the  rate  at  which  monomeric  silicon  collides  with  other  monomers 
or  existing  clusters  to  create  nuclei  of  supercritical  size.  Thus,  the 
description  of  the  nucleation  process  requires  the  determination  of  the 
molecular  silicon  vapor  concentration,  the  supersa tura tion  ratio,  and  the 
nucleation  rate  from  the  existing  theories. 

The  supersaturation  ratio,  SS,  is  defined  by  the  quotient  of  the  actual 
partial  pressure  of  Si(g)  in  the  reactant  gas  stream  and  by  the  equilibrium 
vapor  pressure  of  Si°(g) 

p 

SS  *  .  (2) 

PSi°(g) 

Assuming  that  SiH4<g),  Si(g)  and  H2(g)  are  the  predominant  vapor  species,  the 
equilibrium  concentration  of  Si(g)  in  the  reactant  gas  stream  is  given  by  the 
reaction: 

SiH4(g)  ♦  Si(g)  +  2H2(g)  (3) 

with  an  equilibrium  constant 


2 


Psi(g>  ph2 

kl  “  P 

SiH4 

Pg^0,  is  determined  by  the  reaction: 

(A) 

The  equilibrium  vapor  pressure 

Si(s)  -  Si(g) 

(5) 

with  an  equilibrium  constant: 

k2-  Pglo  • 

(6) 

With  known  equilibrium  constants,  P^^  ^  and  Psj°(g)  can  determined. 

Figure  IV. 5  shows  the  effects  of  temperature  and  the  extent  of  gas  phase 
decomposition  on  the  supersa tura tion  ratio  for  a  reaction  zone  initially 
comprised  of  pure  silane  at  2  x  lO^Pa  <0.2  atm).  For  all  temperatures  the 


57 


NUMBER  DENSITY  x  PATHLENGTH  { 


Figure  IV. 4,  Product  of  the  particle  number  density  and  optical  pathlength 
as  a  function  of  distance  above  the  inlet  nozzle. 


Figure  IV. 5.  The  effect  of  temperature  and  percent  reaction  on  the 
supersaturation  ratio  for  silane  pyrolysis. 


supersaturation  ratio  decreases  sharply  from  infinity  as  the  reaction  begins 
to  proceed. 

★ 

Both  the  nucleation  rate,  J,  and  the  critical  radius  size,  r  ,  in  a 
supersaturated  vapor  can  be  approximated  by  the  application  of  homogeneous 
nucleation  theory.  The  critical  size  nucleus  is  related  to  the 
supersaturation  ratio  by  the  Gibbs-Kelvin  equation:21 
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Si 


kTinSS 


(7) 


where  a  is  the  interfacial  energy  between  the  solid  and  vapor  and  uc  is  the 
molecular  volume  of  silicon.  The  free  energy  of  formation  of  the  critical 
size  nucleus  is  given  by: 
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Using  Volmer-Weber  theory,22  the  equation  describing  the  homogeneous 
nucleation  rate  is: 


(8) 
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where  a  is  the  condensation  coefficient  and  m  is  the  mass  of  a  silicon  atom. 
The  nucleation  rate  is  dependent  on  both  the  thermodynamic  driving  force  for 
the  phase  transformation  and  kinetic  factors  related  to  the  collisions  of 
vapor  phase  silicon  atoms  with  critical  size  silicon  nuclei. 

Table  IV. 2  summarizes  nucleation  rate  calculations  for  a  reactant  gas 
stream  Initially  composed  of  pure  silane  gas  at  2  x  104Pa  (0.2  atm).  The 
critical  size  nucleus  contains  only  a  few  silicon  atoms  as  a  consequence  of 
the  extremely  large  supersaturation,  for  example,  SS  =  7.8  x  1049  at  300K. 
Consequently,  there  is  essentially  no  thermodynamic  barrier  for  nucleation. 
Rather,  the  nucleation  rate  is  primarily  controlled  by  the  collision  rate 
between  silicon  atoms  and  subcritical  nuclei;  at  low  temperatures  the  rate  is 
low  because  the  equilibrium  concentration  of  silicon  vapor  is  extremely 
small. 


The  actual  nucleation  rates  can  be  estimated  from  the  scatter-extinction 
measurement.  The  calculated  particle  number  densities  ranged  from  1.1  x  1011 
to  2.2  x  1012  nuclei/cm3.  The  reactant  gas  traveled  from  the  inlet  nozzle  to 
the  position  where  the  particles  were  first  observed  in  less  than  10”2s  with 
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TABLE  IV. 2 


HU CLEAT ION  CALCULATIONS  FOR  A  REACTION  ZONE  INITIALLY 
COMPRISED  OF  0.2  ATM.  OF  PURE  SILANE 


Temp 

(°K) 

PSi 

(atm) 

P  ** 

Si 

(atm) 

In  SS 

* 

r 

(nm) 

★ 

AG 

j> 

J 

(/// cmJ*s 

300 

1.65xl0_21 

2.12xl0"71 

114.9 

0.075 

2. 15xlo"20 

7.4xl0_1 

400 

2.02xl0"15 

8.77xl0"55 

90.6 

0.072 

1.94xlO"20 

3.5xl0~“ 

500 

9.64xl0~12 

5.16x10"“° 

65.1 

0.080 

2 ,40xl0"2  ° 

7.4xl03 

•600 

2.82xl0_9 

3.59xl0-32 

52.7 

0.082 
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7.6x10® 
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2 . 70xl0_2  0 

2.4xl012 

800 
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1.7xl017 
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28.1 
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7.5xl018 
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1.25xl0"3 
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1 . 8xl02  0 
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1.30xl0_12 
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19.5 

0.103 
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1.6xl022 
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5.55xlO"20 

3.2xl02  3 

all  flow  conditions;  a  time  that  may  be  significantly  longer  than  the  period 
in  which  nucleatlon  actually  occurs.  Therefore,  the  actual  SI  nucleation 
rates  are  In  excess  of  101**  nuclei/cm3 «s. 

It  was  not  possible  to  observe  the  actual  nucleatlon  events  with  the 
scatter-extinction  technique10  since  the  minimum  detectable  particle  diameter 
is  approximately  8  nm  with  a  HeNe  light  source.  Also,  the  temperature  at 
which  nucleatlon  begins  could  not  be  measured  by  optical  pyroraetry  as  a 
consequence  of  the  low  radiance  of  small  particles  (~  1.0  nm)  having  true 
temperatures  below  900°C.  The  lowest  temperature  at  which  particles  were 
first  observed  ranged  from  approximately  900°C  to  1100°C.  It  is  believed 
that  the  nucleatlon  actually  occurred  at  temperatures  significantly  below 
these  measured  temperatures:  probably  between  600  and  900°C. 

The  nucleatlon  rates  calculated  with  the  Volmer-Weber  theory, 

Table  IV. 2,  can  be  compared  to  the  nucleation  rate  estimated  from  the 
scatter-extinction  results.  For  pure  S1H4  at  0.2  atm,  the  calculated 
nucleation  rate  varies  from  1.7  x  1017  nuclei/cm3«s  at  900K  to  2.1  x  1021 
nuclei/cm3«s  at  1200K.  These  values  are  approximately  3  to  7  orders  of 
magnitude  higher  than  the  minimum  nucleatlon  rate  value  estimated  by  scatter- 
extinction  measurements. 

There  are  many  possible  reasons  for  the  discrepancy  between  the  measured 
and  calculated  nucleatlon  rates.  It  was  assumed  that  the  condensation 
coefficient,  a,  was  equal  to  1.0  for  the  calculation  of  the  nucleation  rate 
using  the  Volmer-Weber  theory.  However,  the  condensation  coefficient  is 
always  less  than  unity  for  bulk  samples  and  is  expected  to  be  even  lower  for 
nucleus-size  clusters  because  they  are  unable  to  absorb  the  heat  of 
condensation  to  the  extent  that  the  large  clusters  can.23 

It  was  assumed  that  SiH4(g)  and  Si(g)  were  the  primary  silicon- 
containing  vapor  species  within  the  reaction  zone.  If  other  vapor  species 
(S1H2,  SiHj,  etc)  are  present,  then  the  actual  concentration  of  Sl(g)  will  be 
less  than  that  was  calculated  presuming  the  existence  of  only  SlH4(g)  and 
Sl(g).  This  correction  would  be  quite  large,  because  the  calculated 
nucleatlon  rate  is  dependent  on  the  square  of  the  Sl(g)  monomer 
concentration. 

It  Is  probable  that  the  time  interval  during  which  nucleatlon  actually 
occurs  is  much  shorter  than  the  value  estimated  from  the  gas  velocity  and  the 
onset  of  observed  nucleation.  Also,  the  nucleation  temperature  may  be  lower 


than  that  used  In  the  calculations.  Several  plausible  combinations  of 
nucleation  time  and  temperature,  ranging  from  10_7s  at  1000K  to  10-2s  at 
800K,  yield  calculated  nucleation  rates  that  agree  with  experimental 
observations. 

One  of  the  major  uncertainties  of  the  classical  homogeneous  nucleation 
theory  lies  in  the  assumption  that  clusters  containing  only  a  few  atoms  can 
be  described  in  terms  of  bulk  material  properties.  For  single- component 

systems,  several  workers  have  proposed  that  the  interfacial  energy  increases 
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as  the  cluster  radius  decreases.  *  If  the  nucleating  temperature  Is 
assumed  to  be  1000K,  an  increase  in  the  interfacial  energy  from  0.9  J/m2  to 
1.5  J/m2  would  result  in  approximately  a  10 ^  decrease  in  the  calculated 
nucleation  rate,  bringing  the  calculated  value  into  good  agreement  with  the 
measured  nucleation  rate. 

Finally,  there  are  both  corrections  and  alternatives  to  the  classical 
Volmer-Weber  theory.26  The  actual  steady-state  nucleation  rate  is  always 
lower  than  the  equilibrium  rate  predicted  by  the  Volmer-Weber  treatment.  The 
corrected  steady-state  solution27  generally  introduces  a  nonequilibrium  or 
Zeldovich  factor  into  Equation  9.  For  pure  silicon  at  0.2  atm  and  nucleation 
temperatures  of  900K  and  1200K,  the  Zeldovich  factor  is  calculated  to  be  0.51 
and  0.49  respectively.  Therefore,  this  correction  is  only  partially 
responsible  for  the  discrepancy  between  the  measured  and  calculated 
nucleation  rates.  The  Lothe-Pound  relation28  which  considers  the  statistical 
mechanical  contributions  to  the  free  energy  of  formation  of  the  critical 
nuclei  is  an  important  alternative  to  the  classical  nucleation  theory.  At 
T  *  1000K,  for  a  spherical  cluster  containing  2  to  3  silicon  atoms,  the 
Lothe-Pound  correction  factor  is  approximately  102°.  Thus  the  Lothe-Pound 
theory  is  not  applicable  because  it  predicts  an  even  larger  nucleation  rate. 

In  light  of  these  uncertainties,  it  is  of  little  value  to  pursue  the 
differences  between  measured  and  calculated  nucleation  rates.  The  most 
fundamental  problem  is  one  of  utilizing  bulk  material  properties  to  describe 
clusters  containing  only  a  few  Individual  atoms.  We  accept  the  Volmer-Weber 
model  as  an  adequte  semlquantitatlve  description  of  the  nucleation  process. 

The  number  densities  calculated  from  the  scatter-extinction  measurements 
indicate  that  nucleation  ceases  before  the  reaction  proceeds  to  an 
appreciable  extent.  Murthy  et  al.29  observed  similar  behavior  for  the 
thermal  decomposition  of  silane  in  a  resistance  heated  reaction  tube  furnace. 
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They  found  that  particles  formed  at  the  beginning  of  the  hot  zone  of  the 
reactor  (approximately  700°C)  and  that  very  little  further  nucleation 
occurred  after  the  particles  achieved  macroscopic  size. 

The  rapid  cessation  of  nucleation  is  probably  a  consequence  of  two 
factors.  As  solid  silicon  is  produced,  the  silicon  monomer  concentration 
within  the  reaction  zone  is  reduced;  this  is  evident  from  Figure  IV. 5  which 
shows  that  the  supersaturation  ratio  decreases  sharply  as  the  silane 
decomposition  reaction  proceeds.  Nucleation  effectively  ceases  when  the 
silicon  monomer  concentration  drops  below  a  critical  level.  Also,  once 
nuclei  are  formed,  they  act  as  heterogeneous  sites,  and  growth  by  simple 
condensation  is  favored  over  further  nucleation.8 

IV. 4. B.  Growth 

The  light  scattering  results  show  that  once  nucleation  ceases,  the 
particles  grow  to  their  final  dimensions  via  the  accretion  of  vapor  phase 
molecules;  growth  by  coagulation  was  not  evident  under  the  reaction 
conditions  studied.  The  final  particle  si2e  is  limited  by  the  depletion  of 
the  available  silane,  and  not  by  the  elapsed  growth  time  at  elevated 
temperatures.  Thus,  the  final  size  is  dependent  on  the  number  of  nuclei 
formed  and  the  concentration  of  silicon  atoms  available  to  the  growing  nuclei 
rather  than  growth  rate.  However,  analysis  of  the  growth  process  is 
Important  in  understanding  the  particle  formation  mechanisms. 

The  growth  process  involves  transport  of  a  vapor  species  to  the  surface 
of  the  silicon  nucleus  and  condensation  or  deposition  steps  at  the  surface. 

A  typical  vapor  phase  growth  process  involves  the  following  steps. 

1.  Transport  of  the  reactive  species  to  the  particle  surface. 

2.  Adsorption  of  the  reactants  on  the  particle  surface. 

3.  Formation  of  the  solid  product  either  by  simple  condensation  or  a 
chemical  reaction. 

4.  Desorption  of  any  gaseous  reaction  products. 

5.  Transport  of  the  gaseous  reaction  products  away  from  the  particle 
surface. 

The  overall  growth  rate  will  be  determined  by  the  slowest  step.  Steps  1  and 
5  involve  the  transport  of  material  in  the  vapor,  generally  by  gaseous 
diffusion.  Steps  2,  3  and  4  occur  at  or  on  the  particle  surface;  growth 
limited  by  one  of  these  steps  is  termed  chemical  or  surface  controlled. 


The  rate  laws  for  growth  limited  by  vapor  transport  have  been  derived  by 
various  investigators.30  Ignoring  the  effects  on  the  growth  rate  of  a 
plurality  of  particles  and  the  moving  boundary,  the  volumetric  growth  law  for 
particles  much  larger  than  the  mean  free  path  in  the  vapor  phase  is 


dV  DrpUSi  (PSi(g)~PSi°(g)j 

dt  “  kT 

and  for  particles  smaller  than  the  mean  free  path  is 


dV 

dt 


4ai,rp2u  Sl(PSl(g)~PSi°(g)) 
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where  D  is  the  gaseous  diffusion  coefficient,  r^  is  the  particle  radius,  u 
is  the  molecular  volume,  P 


'Si(g) 


p  •  Si 

is  silicon  partial  pressure  in  the  reactant 


gas  stream,  Is  the  equilibrium  silicon  partial  pressure,  and  m  is 

molecular  mass.  For  all  the  reaction  conditions  studied,  the  mean  free  path 

within  the  reaction  zone  was  greater  than  1  pm.  Because  this  is 

significantly  larger  than  the  particle  size,  it  is  anticipated  that  the 

particle  growth  rate  will  be  represented  by  Equation  11  if  the  growth  rate  is 

limited  by  vapor  transport.  Consequently,  the  radial  growth  rate 

dr/dt  *  (1/4tt  r  2)  dV/dt  will  be  independent  of  the  particle  radius  and 
P 

directly  proportional  to  the  concentration  of  Si(g)  since  the  vapor  pressure 

over  solid  silicon,  P„,0/  \  is  much  smaller  than  the  equilibrium  pressure  of 

oi  (g) 

Si(g)  within  the  vapor  phase.  However,  if  vapor  transport  processes  are 
rapid,  particle  growth  rates  will  be  controlled  by  processes  occurring  at  the 
particle  surfaces. 

Seto3l  studied  the  growth  rate  of  polycrystalline  silicon  from  the 
pyrolysis  of  silane  at  temperatures  where  the  growth  kinetics  are  rate 
limited  at  the  surface.  A  reaction  mechanism  involving  the  chemisorption  of 
silane  onto  the  silicon  surface,  the  decomposition  of  silane  into  silicon  and 
hydrogen,  and  the  desorption  of  hydrogen  explained  the  experimental  results 
satisfactorily.  He  found  that  the  desorption  step  could  be  assumed  to  be 
fast  in  comparison  with  adsorption  and  decomposition  steps  because  the 
hydrogen  gas  produced  by  the  decomposition  reaction  did  not  effect  the 
decomposition  rate. 

Seto  found  that  the  reaction  rate  is  determined  by  the  relative  values 
of  the  Langmuir  adsorption  coefficient  and  the  reaction  rate  constant 


6^4 


where  8, 


kSiH./k-SiH  ’  with  kSlH, 

_  .  4  ,4  .  .  4 


representing  the  silane  absorption 


rate  constant  and  k_slH  being  the  desorption  rate  constant,  and  being  the 
rate  constant  for  the  silane  decomposition  reaction: 

SiH^g)  -  Si(s)  +  2H2(g)  .  (12) 

The  overall  growth  rate31  is  represented  by: 

dv  uSikR6SiH 

Tt  “  1  +  6S1H  [SiHj  (13) 

4 

where  [SiH4J  is  the  molar  concentration  of  SiH^  in  the  gas  stream.  For 
silane  concentrations  ranging  from  5  x  10“7  mole/cm3  to  4  x  10~8  mole/cm3  and 
temperatures  from  650  to  800°C,  Seto  found  that  the  following  expressions  fit 
experimental  observations. 


*  2.0  x  101 3  exp  - 


15200 


kR  =  9.4  x  10"4  exp  - .  (15) 

Equations  10,  11  and  13  predict  specific  functional  dependencies  between 
growth  rate  and  variables  such  as  temperature,  particle  size,  and  reactant 
concentration.  Our  synthesis  experiments  permitted  the  identification  of  the 
rate-controlling  step.  Particle  diameter,  number  density,  and  temperature 
were  measured  directly  by  the  scatter-extinction  technique.  Radial  growth 
rates  were  determined  from  data  like  that  shown  in  Figure  IV. 3.  The 
concentration  of  silicon  remaining  in  the  vapor  phase  was  calculated 
throughout  the  reaction  zone  by  subtracting  the  local  amount  of  silicon 
existing  in  solid  particles  from  the  inlet  concentration. 

Figure  IV. 6  shows  the  calculated  concentration  of  silicon  in  the  vapor 
phase  as  a  function  of  position  in  the  reaction  zone  for  Run  654S,  in  which  a 
2:1  Ar:SiH4  gas  mixture  was  used.  It  is  evident  that  the  reaction  occurred 
in  a  continuous  manner  and  was  essentially  complete  while  the  particles  and 
gases  were  still  in  the  laser  beam.  The  concentration  profiles  for  the  other 
runs  were  similar  to  that  shown  in  Figure  IV. 6. 

Figure  IV. 7  plots  the  natural  logarithm  of  the  growth  rate  versus  the 
natural  logarithm  of  the  concentration  of  silicon  within  the  reaction  zone 
for  Runs  630S  and  654S.  The  growth  rates  were  determined  from  the  slopes  of 
the  curves  In  Figure  IV. 3  assuming  that  the  velocity  remains  constant  in  the 
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Figure  IV. 6.  Calculated  silicon  concentration  in  the  reaction  zone  as  a 
function  of  distance  above  the  inlet  nozzle  for  Run  654S. 


Figure  IV. 7.  The  logarithm  of  the  growth  rate  calculated  from  the 

scattering-extinction  results  versus  the  logarithm  of  the 
silicon  concentration  in  the  reaction  zone.  (a.  Run  630S. 
b.  Run  654S.) 
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spatial  intervals  examined.  This  assumption  introduces  little  error  since 
these  intervals  are  small.  For  both  runs,  the  reaction  temperature  is 
relatively  constant  in  the  regions  plotted.  Therefore,  the  change  in  the 
growth  rate  is  expected  to  be  directly  related  to  the  change  in  the  vapor 
concentration.  The  slopes  of  these  curves  are  in  reasonable  agreement  with 
the  value  (1.0)  predicted  by  Equation  11  for  transport  rate-controlled 
growth. 

This  conclusion  agrees  with  findings  of  others.  At  temperatures  above 
900°C,  Bloem32  determined  that  growth  rate  was  dependent  on  vapor  transport 
for  chemical  vapor  deposition  of  silicon  from  silane;  at  lower  temperatures, 
the  surface  reaction  kinetics  determined  the  overall  growth  behavior. 

Everstyn  et  al.33  also  found  the  deposition  from  silane  at  1050°C  to  be  a 
transport  limited  process. 

It  is  possible  to  estimate  the  condensation  coefficient  by  comparing  the 
growth  rate  measured  from  light  scattering  with  that  predicted  by 
Equation  11.  Table  IV. 3  shows  the  growth  rate  calculated  using  Equation  11 
for  Runs  630S  and  654S  (see  Table  IV. 1  for  reaction  conditions).  The  partial 
pressure  of  SiH4(g)  and  SI(g)  was  determined  from  the  amount  of  silicon  in 
the  vapor  phase  (estimated  from  Figure  IV. 6  and  a  similar  curve  for  Run  63 OS) 
and  equilibrium  thermodynamics.  The  silicon  particles  can  grow  by  two 
mechanisms:  through  gas  phase  decomposition  of  the  silane  according  to 

Equation  1  and  condensation  of  the  Si(g)  molecules  onto  the  growing  particle 
or  by  the  adsorption  of  SiH4(g)  molecules  and  the  decomposition  of  adsorbed 
SiH4  to  form  Si(g)  and  H2(g)  which  must  be  desorbed  from  the  particle 
surface.  Therefore,  the  growth  was  calculated  by  assuming:  (1)  only  the 
adsorption  of  silicon  vapor  leads  to  particle  growth  and  (2)  both  silicon  and 
silane  vapor  molecules  are  adsorbed  during  growth.  In  both  cases  a 
condensation  coefficient  of  unity  was  assumed. 

It  is  evident  that  the  experimentally  observed  growth  rate  Is  2  to  3 
orders  of  magnitude  lower  than  the  calculated  growth  rate.  The  effects  of  a 
plurality  of  particles  and  of  the  moving  growth  surface  result  in  negligible 
corrections  to  the  calculated  growth  rate.8  Therefore,  it  is  assumed  the 
growth  is  limited  by  the  bombardment  of  vapor  molecules,  which  have  a 
condensation  coefficient  significantly  less  than  unity;  the  condensation 
coefficient  required  to  bring  the  experimental  and  calculated  results  into 
agreement  is  shown  in  Table  IV. 3.  If  it  is  assumed  that  both  SlH4(g)  and 


TABLE  IV. 3 


COMPARISON  OF  EXPERIMENTAL  GROWTH  RATE 
AND  CALCULATED  GROWTH  RATE  (TRANSPORT  LIMITED) 


Distance  From  P_.„ 
SiH. 


Nozzle  (mm) 


Run  654S 


(Nt/m2  ) 


Growth  Rate  (nm/sec) 
Measured  Calc. A  Calc.B 


Si  Si+SiH, 


4 . lxlO3  6.7xl02  4.5xl03  1.9xl05  l.lxlO6  0.02  0.004 
1.2xl03  9 . lxlO2  2.4xl03  2.4x10s  5.7x10s  0.01  0.004 
6 . lxlO2  8.0xl02  1 .6x10 3  2.2x10s  3.8x10s  0.008  0.004 
4.9xl02  4 . 8xl02  7 .OxlO2  2.5x10s  2.6x10s  0.005  0.003 
7 .OxlO2  2 .6xl02  4 . 8x1 02  7.4xl04  2.7x10s  0.007  0.002 
4.6xl02  1 .OxlO2  2. OxlO2  2.9x10“  1.4x10s  0.007  0.002 


1 .  lxlO4  2 . 2x10 3  1 . 2X101*  6,1x10s  3,8xl06|  0.02  0.003 
7.4xl03  2 . 2x10 3  8.2xl03  6.1x10s  2.6xl06  0.01  0.003 
5 . 2x10 3  2. lxlO3  4 . 2xl03  5.7x10s  1.9xl06  0.008  0.002 
4.7xl03  1 . 2x1 0 3  3.6xl03  3.4x10s  1.7xl06  0.01  0.002 


A  calculates  growth  assuming  only  Si(g)  atoms  are  absorbed. 

B  calculates  growth  assuming  S i ( g )  and  SiH  atoms  are  absorbed 


Si(g)  species  are  adsorbed  on  the  particle  surface,  the  condensation 
coefficient  is  essentially  constant,  ranging  from  0.002  to  0.004.  For 
silicon  adsorption,  the  calculated  condensation  coefficient  is  larger  and 
exhibits  a  larger  range.  The  possible  existence  of  the  vapor  species  SiH^ 
and  SiH3  has  been  ignored.  If  they  are  present  in  the  vapor  phase  and  also 
contribute  to  the  particle  growth,  the  theoretical  growth  rate  will  be 
approximately  equal  to  that  calculated  above  for  SiH4(g)  and  Si(g)  adsorption 
since  the  total  concentration  of  silicon  containing  species  in  the  vapor 
phase  will  be  equivalent. 

IV. 4. C.  Control  of  Particle  Size,  Shape  and  Size  Distribution 

This  powder  synthesis  process  was  investigated  with  the  goal  of 
achieving  control  of  particle  size,  shape,  and  size  distribution.  It  is 
essential  that  resulting  powders  be  uniform  in  size,  equixial  and  small  in 
diameter. 

Since  the  particle  growth  is  limited  by  vapor  transport  and  not  a 
surface  reaction,  a  thermodynamic  equilibrium  can  be  assumed  to  exist  at  the 
particle  surface.  Therefore,  the  growth  rate  will  not  be  dependent  on  the 
crystallographic  nature  of  the  surfaces  and  particles  will  tend  to  grow 
isotropically,  resulting  in  equiaxed  shapes.34  All  powders  examined  by 
electron  microscopy  were  found  to  be  spherical  in  shape. 

The  early  cessation  of  homogeneous  nucleation  during  the  reaction 
results  In  a  large  number  of  independent  nuclei  that  grow  to  their  final 
dimensions  by  the  depletion  of  the  available  silane  gas.  The  general 
separation  of  nucleation  and  growth  processes  tends  to  produce  uniform  size 
particles.  However,  there  Is  a  short  period  of  time  during  which  nucleation 
and  growth  probably  occur  simultaneously,  which  results  in  an  initial 
distribution  of  particle  sizes.  It  can  be  ascertained  whether  the  difference 
in  the  particle  sizes  becomes  accentuated  or  whether  the  system  tends  toward 
monodlsperslvlty  with  further  growth. 

If  all  the  nuclei  are  exposed  to  the  same  temperature  and  concentration 
profiles,  the  differences  in  the  particle  growth  should  result  only  from  the 
Initial  differences  in  particle  size  when  nucleation  ceased.  When  growth  is 
limited  by  vapor  transport,  Firsch  and  Collins35  showed  that  the  growth  rates 
of  particles  with  different  radii  are  related  by  the  expression 


dr/dt  =  y  +  r1 
dr'/dt  y  +  r 


(16) 


where  r  and  r*  are  the  radii  of  particles  that  had  radii  equal  to  r0  and  ri 
at  the  time  when  nucleatlon  ceased,  and  y  =  p/a  where  p  is  on  the  order  of 
the  mean  free  path  and  a  is  the  probability  that  Impinging  molecules  will  be 
adsorbed.  If  y  »  r,  the  particles  will  grow  at  equal  rates.  Even  if  y  >  r 
throughout  the  entire  growth  period,  the  initial  differences  in  the  particle 
sizes  will  become  insignificant  as  the  particles  approach  their  final  sizes, 
because  the  nucleatlon  process  terminated  early  in  the  growth  process.  If 
with  continued  growth  r  >  y,  the  relative  growth  rate  will  become  inversely 
proportional  to  the  ratio  of  the  radii  tending  to  suppress  the  initial 
differences  in  the  relative  particle  radii.  Therefore,  it  may  be  concluded 
that  a  synthesis  process  characterized  by  early  cessation  of  homogeneous 
nucleatlon  and  transport  limited  growth  will  result  in  a  particle  size 
distribution  that  tends  towards  monodispersity .  For  all  synthesis 
conditions,  the  final  particle  size  distribution  was  always  somewhat 
polydispersed.  This  result  is  probably  a  consequence  of  temperature  and 
compositional  gradients  within  the  reaction  zone. 

The  large  thermodynamic  driving  force  for  the  production  of  solid 
silicon  from  silane  gas  leads  to  the  complete  utilization  of  the  silane. 

Final  particle  size  is  determined  by  the  amount  of  silane  available  for  each 
growing  nuclei.  Assuming  100%  conversion,  a  material  balance  gives  the 
following  relation  for  the  final  particle  size: 

,  C„M  1/3 

d  =  -%)  (16) 
where  d  is  the  average  diameter  of  the  particles  formed,  C0  is  the  Initial 
concentration  of  silane  in  the  reaction  zone,  N  is  the  number  of  nuclei 
formed  per  unit  volume  of  reacting  gas,  M  is  the  molecular  weight  of  silicon, 
and  p  is  the  density  of  the  product.  The  particle  size  can  be  controlled  by 
manipulating  the  concentration  of  the  reactant  gas  and  the  number  of  nuclei 
produced. 
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IV. 5.  CONCLUSIONS 


It  was  found  that  the  laser  synthesis  process  produces  powder  with 
nearly  ideal  characteristics.  The  particle  characteristics  can  be 
systematically  and  reproducibly  varied  by  the  manipulation  of  the  laser 
intensity,  reaction  pressure,  dilution  of  the  reactant  gas,  and  reactant  gas 
flow  rates.  The  conversion  efficiency  to  solid  particles  is  nearly  100%. 

The  silane  decomposition  reaction  is  thermal  in  nature  with  the 
temperature  of  the  reactant  gas  determined  by  the  laser  intensity,  the 
effective  absorption  coefficient  at  the  particular  emission  line,  the  heat 
capacity  of  the  reactant  gas  mixture,  and  the  residence  time  in  the  laser 
beam. 

Gas  phase  nucleation  is  initiated  at  the  beginning  of  the  reaction  zone 
and  ceases  just  after  the  reaction  is  initiated.  The  cessation  is  a  result 
of  the  depletion  of  silicon  vapor  reducing  the  supersaturation  as  well  as  the 
kinetic  favoring  of  particle  growth  by  the  accretion  of  vapor  molecules  over 
further  nucleation.  The  observed  nucleation  rate,  >1014  nuclei/cm3«sec,  is 
in  reasonable  agreement  with  the  calculated  rate  using  Volmer-Weber  theory. 

The  growth  was  determined  to  be  rate  limited  by  the  transport  of  vapor 
species  to  the  particle  surface.  The  experimentally  observed  growth  rate  is 
2  to  3  orders  of  magnitude  lower  than  the  growth  rate  calculated  from  the 
vapor  transport  model  assuming  a  condensation  coefficient  of  unity.  A 
condensation  coefficient  of  0.002  to  0.02  brings  the  theoretical  growth  rate 
into  agreement  with  the  observed  growth  rate.  The  nature  of  the  growth 
process  results  in  spherical  particles  with  narrow  distribution  of  particle 
sizes. 

Under  the  synthesis  conditions  investigated,  particle  agglomeration  did 
not  occur  within  the  reaction  zone.  Chainlike  agglomerates  formed  as  the 
particles  travelled  to  the  collection  filter  or  as  they  were  collected. 
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CHAPTER  V 


DISPERSION  OF  SILICON  IN  N-PROPANOL 

by 

G.  Garvey  and  J.  S.  Haggerty 

ABSTRACT 

Dispersions  of  laser-formed  silicon  powders  in  n-propanol  were 
evaluated.  Parameters  studied  included  boron  doping,  fraction  solids  and 
water  concentration  in  the  n-propanol.  FTIP.  was  used  to  evaluate  the 
chemical  composition  of  the  powder.  Powder  agglomeration  was  evaluated  by 
TEM  and  photon  correlation  spectroscopy;  macroscopic  dispersion  stability  was 
measured  by  sedimentation  rates.  Microelectrophoretic  mobility  was  used  to 
evaluate  surface  potentials. 

The  results  are  interpreted  in  terms  of  the  DLVO  model.  In  general, 
coulombic  repulsion  does  not  provide  an  adequate  potential  barrier  because 
the  double  layer  is  very  large  in  the  alcohol.  Water  has  an  adverse  effect 
on  stability.  It  may  be  possible  to  achieve  stable  dispersions  with  ionic 
additions;  but,  steric  stabilization  by  adsorbed  macromolecules  appears  more 
likely  to  succeed. 

V.l.  INTRODUCTION 

The  objective  of  the  dispersion  task  of  this  program  is  to  find  means  of 
dispersing  nonpolar,  covalent  particles  in  a  manner  that  will  permit  defect 
free,  high  coordination  number  green  bodies  to  be  formed.  Dispersion  must  be 
accomplished  in  liquids  and  with  additives  that  will  not  degrade  the 
properties  of  densified  bodies.  Superior  microstructural  characteristics  and 
reduced  sintering  cycles  have  been  demonstrated  with  TiOjl  with  green  bodies 
made  from  the  special  dispersions  that  are  possible  only  with  uniform 
diameter  particles.  This  research  strives  to  achieve  similar  results  with 
Si,  SIjN4  and  SiC  powders  produced  by  the  laser  heated  gas  phase  synthesis 
process. 

It  was  postulated2  and  now  has  been  shown3  that  essentially  defect  free 
green  ceramic  bodies  can  be  formed  from  ordered  dispersions  of  ceramic 
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particles.  The  particles  must  conform  to  a  rigid  set  of  criteria  for  ordered 
dispersions  to  be  feasible;  principally  they  must  be  uniform  in  diameter  and 
absolutely  free  of  agglomerates.  Conformance  to  other  criteria  requires 
specific  particle  sizes,  purities  and  stoichiometries.  The  laser  heated  gas 
phase  synthesis  process  generally  appears  capable  of  producing  powders  with 
these  characteristics. 

Establishment  of  ordered  dispersions  with  minimum  interparticle 
distances  requires  precise  manipulation  of  interparticle  forces  thereby 
maximizing  the  volume  fraction  of  solids.  The  theoretical  understanding  of 
these  interparticle  forces  in  nonaqueous  liquids  is  far  less  developed  than 
for  ionic  particles  in  aqueous  liquids.  Thus,  our  research  has  begun  with 
empirical  screenings  and  fundamental  property  measurements  to  permit 
definition  of  guidelines  for  feasible  approaches. 

The  research  supported  by  this  program  has  focused  on  Si  powders  in 

alcohols,  building  on  the  research  carried  out  in  programs  previously 

\5 

supported  by  ONR  and  ARO.  Dispersion  research  with  laser  synthesized  SiC 

powders  is  being  supported  by  DOE6  and  NASA7.  Combined  they  will  extend 
MIT's  generic  research  in  presintering  science  to  this  important  class  of 
ceramics. 

Silicon  was  selected  for  these  studies  for  several  important  reasons. 

Its  higher  overall  purity,  our  ability  to  control  particle  size  over  wider 
range,  the  absence  of  crystalline  polytypes  and  the  nonexistence  of  the 
stoichiometry  issue  make  Si  a  simpler  system  than  the  other  two  materials. 

It  is  an  excellent  choice  for  modelling  dispersion  of  covalent  materials  in 
nonaqueous  liquids.  Also,  the  material  has  great  technical  importance  for 
reaction  bonded  silicon  nitride  and  reaction  bonded  silicon  carbide. 

Beyond  the  obvious  adverse  effect  flaws  have  on  strength,  highly  perfect 
green  microstructures  are  particularly  important  for  RBSN  products.  Specific 
surface  areas  and  interparticle  neck  sizes  control  the  uniformity  and  rates 
of  nitridation;  both  issues  have  been  troublesome  in  RBSN  processing.  In 
principle  it  should  be  possible  to  make  nearly  fully  dense  RBSN  bodies 
without  any  dimensional  change  from  the  ordered,  green  Si  bodies.  Also,  the 
absence  of  large  diameter  interparticle  voids  will  improve  the  oxidation 
properties  of  RBSN®,  permitting  superior  high  temperature  properties 
exhibited  in  nonoxidizing  atmospheres  to  be  retained  in  oxidizing 
a  tmospheres. 


V.2.  THEORETICAL  CONSIDERATIONS 


Particles  can  exist  in  stable  suspension  if  attractive  and  repulsive 
forces  are  equal  in  magnitude  at  a  noncontacting  distance  and  repulsive 
forces  exceed  attractive  forces  for  interparticle  distances  less  than  the 
equilibrium  distance.  Thermally  induced  oscillations  cause  interparticle 
positions  to  change  continuously  about  the  equilibrium  position;  instability 
results  if  the  repulsive  force  barriers  are  exceeded.  Van  der  Waals  forces 
attract  particles,  steric  barriers  repel,  and  coulombic  forces  can  either 
attract  or  repel.  Typically  coulombic  forces  are  used  to  repel  particles 
against  van  der  Waals  attraction.  Each  force  exhibits  a  different  functional 
dependence  on  interparticle  spacing  and  all  are  extremely  sensitive  to  the 
properties  of  the  specific  materials  (powders,  solvents,  solutes...). 

The  particle  interactions  are  generally  described  in  terms  of  an  energy 
which  is  the  integral  of  the  force-distance  relationship  corresponding  to 
bringing  two  particles  from  an  infinite  separation  to  a  specific 
interparticle  spacing.  The  total  interaction  energy  is  the  sum  of  the 
individual  contributions.  The  relative  stability  of  a  dispersion  is 
dependent  on  the  barrier  height- to- thermal  energy  ratio. 

A  charged  particle  in  a  partially  ionized  medium  will  attract  Ions  of 
the  opposite  sign  which  can  either  adsorb  (Stem  layer)  onto  the  particle 
surface  or  remain  in  the  medium  surrounding  the  particle.  These 
redistributed  ions  constitute  the  electrical  double  layer.  When  two  such 
particles  approach  each  other,  the  diffuse  clouds  of  counter  ions  interact  by 
repelling  each  other.  The  structure  of  the  double  layer  is  governed  by  the 
type  and  concentration  of  ions  in  solution  (in^),  and  the  dielectric  constant 
of  the  medium  (e)  and  the  sign  and  magnitude  of  the  surface  charge.  The 
double  layer  thickness,  1/k,  Is  a  measure  of  the  structure  of  the  double 


layer.  It  is  both  the  thickness  of  an  equivalent  capacitor  and  the  distance 
at  which  the  electrical  potential  in  the  medium  drops  to  1/e  of  the 
particle's  surface  potential  (v0).  This  thickness  is  given  by:9 


l/<  -  (  - - )1/2 

4*  eZ&i^z^2 


where  e  is  the  electronic  charge,  z^  is  the  valence  of  the  ion  i,  and  kT  has 
its  normal  meaning. 
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The  ionic  concentration  in  a  pure  aqueous  medium  is  approximately 
10”  7m  at  room  temperature  (defined  by  its  self  dissociation  constant). 
Electrolytes  are  normally  added  to  dispersions  to  control  surface  charge; 
typical  ionic  concentrations  range  from  10“5  to  10-1m.  These  values  result 
in  double-layer  thicknesses  of  0.1  pm  and  0.001  pm  respectively.  In  pure  n- 
propanol,  the  ionic  concentration  is  approximately  10”9m;  the  resulting 
double- layer  thickness  is  approximately  2  pm  as  shown  in  Figure  V.l.  The 
addition  of  an  organic  salt  which  dissociates  to  a  high  degree  in  n-propanol 
may  be  useful  in  condensing  the  double  layer  thickness  (l/<),  although, 
surface  charge  neutralization  may  occur. 

Although  the  equations  that  govern  the  coulombic  interactional  energies 
have  not  been  solved  explicitly,  solutions  exist  for  the  limiting  cases  of 
large  double  layers  (<a  <  2)  and  small  double  layers  (xa  >  5)  where  a  = 
particle  radius.  Organic  dispersions  fall  into  the  former  case  with  <a  = 
0.01  for  0.1  pm  diameter  Si  particles  in  n-propanol.  For  this  case,  the 
coulombic  repulsive  potential  (V^)  as  a  function  of  interparticle  distance 
(a)  is  given  by:9 


ea 


2*2 


V  = 


Be- 


•c  A 


R  A  +2a 

where  B  is  a  geometrical  correction  factor  which  ranges  from  0.6  to  1.0.  An 
error  of  less  than  one  percent  results  from  ignoring  the  exponential  term  for 
<A  <0.01  (i.e.  A  <  500A).  The  denominator  is  dominated  by  particle  radius 
for  close  interparticle  spacings,  so  the  potential  barrier  does  not  rise 
rapidly  as  the  particles  approach. 

There  are  three  mechanisms  for  generating  surface  charge  and  thus 
surface  potential,  in  the  systems  under  consideration.  These  are: 

•  Distribution  of  electrons  within  the  semiconducting  particles: 

Surface  states  can  adsorb  electrons  or  holes  resulting  in  the 
development  of  a  surface  charge. 

•  Adsorption  of  dipolar  molecules:  The  adsorption  of  neutral  dipolar 
molecules  from  solution  can  effect  the  space  charge11  distribution  in 
the  double  layer. 

•  Surface  reaction:  Any  reaction  that  results  in  the  production  and 
desorption  of  ions  at  the  surface  will  generate  a  surface 
charge;11”15  e.g.,  Bronsted  and  Lewis  mechanisms. 
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Surface  charge  in  silicon  is  expected  from  some  combination  of  all  three 
of  these  mechanisms.  The  magnitudes  of  each  are  interdependent.  Electronic 
distribution  can  be  most  easily  manipulated  by  varying  the  dopant  level  in 
the  silicon  and  by  the  specific  electrical  characteristics  of  the  solvent. 
Dipolar  effects  scale  with  the  dipole  moment  of  the  solvent.  Surface 
reaction  is  manipulated  by  controlling  the  solid  surface  and  solvent 
chemistries.  The  sign  of  each  contribution  can  be  either  positive  or 
negative.  Because  the  potential  decays  slowly  from  the  surface  in  nonaqueous 
liquids,  the  surface  potential  can  be  equated  to  the  ^-potential. 14  The 
magnitude  and  sign  of  the  net  surface  potential  can  be  measured  by 
electrophoretic  experiments . 1 1 

The  van  der  Waals  potential  energy  between  two  equal  diameter  particles 
of  the  same  material  is  given  by: 

V--  — 

12  x2  +  2x 

where  x  is  A/2a  and  A  is  the  Hamaker  constant.  Because  the  denominator  goes 
to  zero  as  A  goes  to  zero,  these  forces  become  large  for  small  interparticle 
spacing.  Literature  values  of  A  used  in  the  following  calculations  are 
summarized  in  Table  V. 1. 

Steric  interference  causes  repulsive  forces  between  interacting 
particles  with  adsorbed  macromolecular  layers.  The  magnitude  of  the 
repulsive  force  depends  on  the  macromolecule  size,  adsorption  density  and 
conformation.  Most  of  the  theories  for  steric  stabilization  can  be 
classified  in  either  entropic  or  osmotic  categories.  Decreased 
configurational  entropy  caused  by  compressed,  interpenetrating  adsorbed 
molecules  is  responsible  for  repulsive  forces  in  the  entropic  theory. 15-2C 
For  the  osmotic  theory,  the  repulsive  forces  result  from  a  combination  of 
positive  enthalpic  change  and  reduced  configurational  entropy  during 
Interparticle  collisions.21”24  Entropic  and  osmotic  repulsion  theories  have 
been  tested  by  many  authors  and  both  have  been  found  to  apply  for  specific 
systems. 

Although  there  is  relatively  little  guidance  for  steric  stabilization  in 
the  Si/n-propanol  system,  this  approach  must  be  considered  very  important  in 
light  of  the  diffuse  nature  of  coulombic  repulsive  potentials.  Steric 
repulsive  forces  operate  at  interparticle  distances  approximately  the  same  as 
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TABLE  V.l 


HAMAKER  CONSTANTS 


Material  Interaction  A  (J  x  IQ20) 

Sl-Sl  25.6* 

n-Propanol-n-Propanol  5.0+ 

H20-H20  4.4* 

*  J.  Visser,  "On  Hamaker  Constants",  Adv.  Coll.  Inter.  Scl.,  3,331-63,  1972 

+  B.  Vincent,  "The  van  der  Waals  Attraction  Between  Colloidal  Particles 
Having  Adsorbed  Layers.  II",  J.  Coll.  Inter.  Scl.,  42,  270-85,  1973. 


LOS  IONIC  CONCENTRATION  C»> 


Figure  V.l.  Double  layer  thickness  as  a  function  of  Ionic  concentration, 
c  -  20.4;  T  -  300 °K . 


I  — . . . ... .. _ _ 
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Figure  V.2.  Schematic  of  interparticle  energies  for  a  typical  stable  aqueous 
dispersion. 


1  I*  IM  IMS 


MXKIMUn  DISTANCE  CA> 

Figure  V. 3.  Interparticle  energies  for  charged  silicon  particles  in 

n-propanol  as  a  function  of  interparticle  distance;  particle 
radius  ■  4004;  surface  charge  ■  80  mV;  [ion]  “lx  10“9  m  at 


Figure  V.4.  Interparticle  energies  for  silicon  in  n-propanol  as  functions  of 
interparticle  distance  and  surface  charge.  Radius  =  400A ; 

[ion]  =  1  x  10"9  m;  T  =  300°K. 


Figure  V. 5.  Interparticle  energies  for  silicon  in  n-propanol  as  functions  of 
interparticle  distance  A  and  particle  radius,  80  mV; 

[ion]  -  1  x  10"9  m;  T  -  300°K. 


the  macromolecule  sizes.  Thus,  they  provide  repulsive  forces  that  operate  on 
the  same  scale  as  the  van  der  Waals  forces.  If  balanced,  their  use  can 
produce  a  stable  dispersion  with  very  small  interparticle  spacings;  a 
critical  feature  for  achieving  a  high  percentage  of  solids  in  green  bodies 
with  small  particles. 

The  van  der  Waals  and  coulombic  potential  terms  included  in  Figures  V.2 
and  V.3  compare  these  interactional  energies  in  a  typical  stable  aqueous 
dispersion  and  in  a  Si/n-propanol  dispersion  respectively.  A  substantial 
potential  barrier  is  developed  in  the  aqueous  example  (Figure  V.2).  If  this 
barrier  is  greater  than  10  kT,  the  dispersion  will  be  stable  to  flocculation 
for  a  period  of  days.  In  contrast,  the  maximum  potential  is  less  than  10  kT 
and  the  potential  barrier  is  very  diffuse  for  the  nonaqueous  example, 

Figure  V.3.  Substantial  changes  in  the  interparticle  distance  result  in 
relatively  minor  changes  in  potential.  For  example,  the  increase  in 
potential  experienced  by  two  particles  as  they  approach  a  separation  of  1001 
from  an  initial  distance  of  1000A  is  small,  on  the  order  of  1  kT.  This 
approach  can  be  driven  by  thermal  energy.  Once  the  particles  reach  a 
distance  of  100,1,  agglomeration  is  energetically  favored.  Such  a  dispersion 
is  not  stable. 

Figures  V.4  and  V.5  show  the  effects  of  surface  potential  and  particle 
diameter  for  Si  particles  in  n-propanol.  These  plots  are  similar  to  those  of 
Figure  V.3;  the  potential  barriers  are  small  and  diffuse.  Varying  the 
particle  size  or  surface  potential  within  these  limits  (0  to  80  mV,  50  to 
800A)  does  not  result  In  stable  dispersions. 

The  Void25  approach  can  be  used  to  predict  the  effects  of  an  adsorbed 
water  layer  on  the  van  der  Waals  potential.  A  brief  discussion  of  this  model 
is  given  by  Parfitt11.  In  short,  the  model  assumes  that  the  total 
van  der  Waals'  potential  of  two  particles  with  adsorbed  layers  is  the  sum  of 
the  interactions  of  each  of  the  phases  with  every  other  phase  present  in  the 
system.  The  adsorbed  layer  is  treated  as  a  sphere  of  liquid  with  a  sphere  of 
the  particle's  geometry  subtracted  from  it.  Figure  V.6  shows  that  for  all 
lnterpartlcle  separations,  adsorbed  water  increases  the  van  der  Waals 
attractive  potential.  The  presence  of  a  20A  layer  of  pure  water  may  change 
the  surface  potential  substantially,  but  is  not  expected  to  affect  the 
diffuse  double  layer  structure  significant]  .  On  the  basis  of  these  results 
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it  is  expected  that  dispersions  of  silicon  in  n-propanol  containing  water 
should  be  less  stable  than  those  without  water. 


In  summary,  interparticle  forces  have  the  same  origins  in  aqueous  and 
nonaqueous  liquids;  however,  their  character  differs.  The  nature  of  the 
coulombic  repulsive  potential  between  silicon  particles  in  pure  n-propanol  is 
the  most  important  difference.  In  contrast  to  aqueous  dispersions,  the  weak, 
diffuse  coulombic  repulsion  does  not  provide  a  sufficient  potential  barrier 
to  prevent  coagulation.  Although  the  approach  is  unproven,  it  may  be 
possible  to  compress  the  coulombic  barrier  (make  less  diffuse)  and  increase 
the  surface  potential  with  easily  dissociated  electrolytes  as  is  done  with 
aqueous  dispersions.  Steric  stabilization  appears  particularly  important  for 
these  dispersions  because  these  forces  act  on  very  small  interparticle 
distances  providing  the  opportunity  to  achieve  stability  with  high  percent 
solids  of  small  particles.  Water  impurities  can  be  anticipated  to  reduce 
stability  through  increased  van  der  Waals  attraction.  Experiments  were 
designed  to  examine  these  issues. 


Figure  V. 6.  Interparticle  attraction  of  silicon  in  n-propanol  with  and 

without  adsorbed  water.  Radius  *  100A;  water  layer  thickness 
-  201;  T  -  300°K. 


V.3.  PROCEDURES 


V.3.A.  IR  Spectroscopy 

The  infra-red  spectra  (IR)  of  two  silicon  powders  (226S  and  227SB)  were 
obtained  using  a  Nicolet  7000  FTIR.  Pressed  pellets  were  prepared  by  milling 
0.1  mg  of  silicon  powder  with  0.4  g  IR-grade  KBr  obtained  from  Fisher 
Scientific.  Pellets  were  pressed  in  a  round  1.27  cm  diameter  Perkin-Elmer 
stainless  steel  die  to  1.4  x  lO^MPa.  A  blank  sample  was  similarly  prepared 
at  the  same  time  for  background  determination.  After  the  spectra  were 
obtained  and  stored,  the  base  line  corrections  were  made  by  subtracting  the 
blank's  spectrum. 

V.3.B.  Specific  Surface  Area,  TEM  and  Photon  Correlation  Determinations 
of  Particle  and  Agglomerate  Sizes 

Specific  surface  areas  were  determined  by  single-point  BET  analysis  with 
a  30%  N2~70%  He  atmosphere.  Equivalent  spherical  diameters  were  calculated 
assuming  a  density  of  2.33  (g/cm3). 

Particles  and  agglomerate  sizes  were  measured  after  powders  had  been 
dispersed.  The  dispersions  were  formed  in  the  following  manner.  A  10  mg 
sample  of  dry  powder  was  weighed  Into  a  Pyrex  scintillation  vial  in  a  dry, 
deoxygenated,  nitrogen-filled  glove  box.  The  vial  was  stoppered  with  a 
rubber  septum  and  removed  from  the  box.  Stock  n-propanol  (10  ml)  was 
injected  Into  the  vial  via  a  polypropylene  syringe.  The  water  content  in  the 
alcohol  was  determined  to  be  ~  250  ppm  by  weight  using  Karl  Fisher  Reagent 
(KFR)  titration.  The  vial  was  shaken  and  sonicated  using  an  Ultrasonics  and 
Heat  Systems  model  W-220  ul trasonica tor  at  a  setting  of  6  for  three  minutes. 
The  exterior  of  the  vial  was  maintained  at  room  temperature  using  a  water 
ba  th. 

For  photon  correlation  (PC)  measurements,  the  dispersion  was  transferred 
to  a  PC  cuvette  and  evaluated  in  a  Malvern  K7025  Photon  Correlator.  Both  the 
transfer  of  the  dispersed  particles  and  the  PC  characterization  were  done 
with  air  exposures.  Measurements  were  Initiated  2  minutes  after  the 
cessation  of  sonication.  For  TEM  measurements,  a  carbon  TEM  film  grid  was 
dipped  Into  the  dispersion  2  minutes  after  sonication  then  dipped  into  pure 
n-propanol.  The  second  dip  prevented  a  nonrepresentative  sample  of  silicon 
dispersion  from  concentrating  on  the  grid  as  the  alcohol  evaporated.  It  is 


nnz>**o*m> 


i 


Figure  V.7. 


FTIR  absorbance  plot  of  sample  226S; 


KBr  pellet. 


Figure  V. 8.  FTIR  absorbance  plot  of  sample  227SB;  KBr  pellet. 
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believed  that  the  silicon  agglomerates  that  remained  on  the  film  were 
representative  of  the  agglomerates  in  the  10  ppm  dispersion  although  this  is 
difficult  to  demonstrate  unequivocably. 

The  effects  of  water  on  the  agglomerate  size  were  determined  by  photon 
correlation.  Three  dispersions  were  formed  as  described  above,  starting  with 
stock  n-propanol  containing  approximately  250  ppm  by  weight  water.  A  small 
amount  of  water  was  added  to  two  of  these  vials  gravimetrically  bringing  the 
water  contents  to  2740  and  13,700  ppm  by  weight.  The  dispersions  were 
ultrasonicated,  transferred  to  PC  vials  and  assessed  for  agglomerate  size  as 
a  function  of  time. 

V.3.C.  Electrophoretic  Mobility  Determinations 

Microelectrophoretic  mobility  was  determined  by  the  method  described  by 
R.  J.  Hunter26  using  the  Rank  Brothers  apparatus.  A  cylindrical  vitreous 
silica  electrophoretic  cell  with  2  mm  ID  was  used.  Applied  fields  ranged 
from  3  to  7  V/cm.  ^-potentials  were  calculated  as  indicated  in  the  cited 
reference.  Typical  values  ranged  from  -30  to  -70  mV.  A  minimum  of  forty 
measurements  were  made  at  two  different  voltages  with  direction  alternated. 
N-propanol  used  in  this  study  was  distilled  after  drying  over  4x  molecular 
sieves  for  3  days  or  was  used  in  the  as-received  condition.  Water  content 
was  determined  by  Karl  Fisher  reagent  titration.  All  glassware  was  dried 
with  a  heat  gun  in  a  N2  glove  box;  the  02  and  H20  content  in  the  box  was 
minimized  by  a  copper  getter.  The  microelectrophoretic  cell  was  loaded  and 
capped  in  the  glove  box. 

A  0.04  wt.%  dispersion  of  silicon  was  prepared  with  water  assayed  n- 
propanol  as  described  above.  This  was  again  diluted  by  a  factor  of  100:1 
with  the  appropriate  n-propanol  solution  to  produce  a  4  ppm  by  weight 
dispersion  which  was  sonicated  as  previously  described.  The  diluted 
dispersion  was  used  for  microelectrophoretic  mobility  determination. 

V.4.  RESULTS  AND  DISCUSSION 

V.4.A.  IR  Study 

The  FTIR  spectrographic  analyses  of  undoped  and  boron  doped  Si  powders 
are  shown  in  Figures  V.7  and  V.8  respectively.  The  spectrum  of  the  boron 
doped  powder  (Figure  V.8,  227-SB)  has  a  very  broad  spectral  roll  centered  at 
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2200  cm-1*.  V.  S.  Vavilow,  et  al.28  have  reported  a  similar,  less  intense 
spectral  feature  for  boron  doped  single  crystal  silicon. 

Two  large  absorption  bands  at  3435  cm-*  and  1645  cm-^  are  apparent  in 
the  spectrum  of  227-SB  that  are  not  present  in  the  spectrum  of  226-S.  0-H 

stretching  vibrations  of  the  water  molecule  lead  to  a  very  broad  and  intense 
band  near  3400  cm-l  and  H-O-H  bending  vibrations  of  molecular  water  lead  to  a 
band  near  1650  cm-!.29  In  contrast,  surface  hydroxyl  in  the  form  of  silanol, 
SiOH,  is  characterized  by  a  very  narrow  band29  near  3750  cm-*.  This  species 
does  not  produce  a  band  near  1650  cm-l.  As  water  adsorbs  on  a  silanol 
surface  the  band  at  3750  cm-1  decreases  in  intensity  and  is  displaced  to 
lower  frequencies29  until  the  peak  is  obscured  by  the  water  spectrum.  The 
observed  results  may  be  explained  by  concluding  that  the  boron  doped  silicon 
is  hydroxylated  and  covered  by  an  adsorbed  layer  of  molecular  water.  Pure 
silicon  has  neither  adsorbed  water  nor  a  hydroxylated  surface  because  neither 
of  the  characteristic  peaks  are  observed. 

Both  spectra  contain  broad,  small  absorbance  peaks  at  the  low  wave 

numbers.  A  peak  centered  at  1100  cm-l  is  most  easily  discerned  in  the  226-S 

spectrum.  Many  absorption  bands  in  this  region  (e.g.  1106  cm-1*,  1136  cm-*, 

1127  cm-1* ,  1130  cm-1  and  1111  cm-*)  have  been  examined  and  attributed  to 

30  31  32 

lattice  oxygen  impurities.  ’  ’ 

Peaks  observed  in  the  420-640  cm-*  region  may  be  due  to  a  Si-O-Si 
symmetric  stretch.33  Also,  Si-H  species  (e.g.  SiH  and  SiH2)  exhibit 
vibration  modes  in  this  region.  Such  species  can  be  anticipated,  being 
products  of  silane  decomposition.  The  appearance  of  a  peak  at  2350  cm-1 
reinforces  the  case  for  the  presence  of  SiH  which  has  a  peak  assigned  in  this 
region. 34 

In  summary,  it  appears  that  both  powders  contain  some  subsurface, 
lattice  oxygen  contamination.  This  implies  that  a  small  amount  of  oxygen 
enters  the  system  during  the  synthesis  process.  Possible  sources  include 
contamination  of  the  reactant  and  shielding  gases  and  leaks.  In  addition  the 
boron-doped  silicon  (227-SB)  has  adsorbed  molecular  water  and  perhaps  some 
surface  silanol  unlike  the  pure  silicon.  This  water  is  probably  introduced 
during  transfer  of  the  powder  after  forming.  A  hydroxylated  surface  is 
expected  to  exist  for  this  powder  but  its  broad  characteristic  peak  is 
obscured  by  water  of  hydration  which  has  a  large  similar  peak  in  this  region. 
The  disposition  of  the  B-doped  silicon  to  adsorb  water  corresponds  to  other 
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Figure  V.ll.  TEM  particle  size  distribution  of  sample  806SB. 
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Figure  V.12.  BET  equivalent  spherical  diameters  versus  TEM  arithmetic  means 
for  eight  silicon  powders. 


independent  observations.  DahlAn35  showed  that  boron  doped  silicon  is  wetted 
by  water  while  nonoxidized,  pure  silicon  powder  is  not  wetted  by  water  and 
floats  on  the  surface.  From  peaks  appearing  at  2350  cm-1,  495  cm"!  and 
570  cm-1,  it  is  concluded  that  SiH  is  also  present  in  the  powders. 

V.4.B.  Particle  and  Agglomerate  Size  Distributions 

Samples  were  prepared  for  particle  size  analysis  as  described  in 
Section  V.3.B.  Particle  size  measurements  were  made  from  TEM 
photomicrographs.  A  minimum  of  400  particles  were  sampled  randomly  from  each 
powder  sample.  The  results  are  reported  in  Figures  V.9,  V.10  and  V.ll. 

Figures  V.10  and  V.ll  reflect  the  observance  of  three  different  types  of 
particles.  Type  one  particles  are  small  (<  300A),  have  a  narrow  particle 
size  distribution  and  commonly  occur  in  agglomerates  of  more  than  twenty.  It 
is  not  evident  from  the  micrographs  whether  these  particlate  assemblages  are 
hard  aggregates  (i.e.  have  sintered  necks).  However  it  is  apparent  that  they 
either  are  not  dispersed  by  the  applied  techniques  or  agglomerate  during  TEM 
sample  preparation.  The  type  two  particles  referred  to  in  Figures  V.10  and 
V.ll,  with  an  average  diameter  of  500A,  commonly  occur  In  chainlike 
agglomerates  some  of  which  appear  to  be  aggregated.  The  number  of  particles 
per  agglomerate,  ~  5,  is  less  than  the  number  for  agglomerates  of  type  one 
particles.  The  type  three  particles  are  large  and  spherical,  occuring 
normally  as  single,  non-agglomerated  particles.  The  weight  fractions  of  each 
type  depend  on  the  individual  powder  sample. 

The  results  of  TEM  particle  size  distributions  were  compared  with  BET 
equivalent  spherical  diameter  measurements.  The  TEM  arithmetic  means  are 
compared  with  BET  equivalent  spherical  diameters  in  Figure  V.12.  As 
expected,  the  TEM  arithmetic  means  are  consistently  lower  than  BET  equivalent 
spherical  diameters  which  weight  larger  particles  more.  Figure  V.13  applies 
the  BET  weighting  factor  to  the  TEM  particle  size  distributions.  There  Is 
very  good  agreement  between  the  two  methods;  however,  the  TEM  values  are 
again  consistently  smaller.  This  small  difference  may  result  from  neck 
formation  between  the  particles  because  simple  particle-particle  contact  does 
not  account  for  all  of  the  discrepancy. 

The  number  of  particles  per  agglomerate  are  presented  in  Figures  V.14 
and  V.15.  Many  of  the  particle  assemblages  are  singlets  or  doublets.  This 
plot  represents  nearly  the  maximum  dispersion  possible  by  the  sample 
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Figure  V.13.  BET  equivalent  spherical  diameters  versus  TEM  weighted  mean 
diameters  for  eight  silicon  powders. 
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Figure  V.14.  TEM  agglomerate  size  distribution  for  sample  226S. 
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Figure  V.17.  TEM  agglomerate  size  distribution  for  sample  226S;  10  ppm 
silicon  in  n-propanol. 
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Figure  V.18.  TEM  agglomerate  size  distribution  for  sample  227SB;  10  ppm 
silicon  in  n-propanol. 
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Figure  V.21.  Log-normal  agglomerate  size  plot  of  data  in  Figure  V.18  sample 
227SB. 
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Figure  V.22.  TEN  particle  size  distribution  for  sample  B004S  showing  three 
types  of  particles.  Bar  length  *  10,0001. 


preparation  procedure,  as  samples  were  taken  only  a  few  minutes  after 
sonication. 

Figures  V.16,  V.17  and  V.18  report  the  results  of  the  TEM  agglomerate 
size  distribution  study.  The  type  of  distribution  must  be  known  to  properly 
characterize  these  results  with  respect  to  statistical  means  and  standard 
deviations.  The  results  were  tested  for  their  agreement  with  log-normal, 
normal,  Weibull,  exponential  and  logistic  distributions.  The  log-normal 
distribution  gave  good  straight  line  agreement  as  indicated  in  Figures  V.19, 
V.20  and  V.21.  Values  of  the  log-normal  means  and  standard  deviations  are 
tabulated  in  Table  V.2.  The  values  of  the  means  and  standard  deviations 
obtained  by  photon  correlation  are  also  tabulated  in  Table  V.2.  Very  good 
agreement  is  found  to  exist  between  the  TEM  and  the  photon  correlation 
distribution  means.  The  standard  deviations  differ  somewhat  more  but  are 
within  acceptable  limits. 

The  effect  of  centrifugal  separation  on  particle  size  and  type  was 
studied  with  powder  B-004-S.  This  powder  was  chosen  because  it  had  a  large 
fraction  of  type  three  particles,  Figure  V.22.  A  0. 1  wt.%  dispersion  in  n- 
propanol  was  prepaced  as  described  in  V.3.B.  The  dispersion  was  centrifuged 
for  a  time  and  speed  calculated  to  settle  2000A  particles  from  the  top  of  the 
container.  The  supernate  was  decanted  and  the  sediment  redispersed  in  fresh 
n-propanol  a  total  of  three  times.  A  TEM  micrograph  of  the  final  product  is 
presented  in  Figure  V.23.  The  majority  of  the  particles  are  large,  spherical 
and  unaggregated.  Also,  there  is  a  narrow  particle  size  distribution  which 
is  necessary  for  the  formation  of  ordered  compacts. 

These  characterization  results  provide  a  self  consistant  description  of 
the  powders.  The  standard  deviation  of  particle  sizes  is  fairly  small, 
usually  less  than  half  the  mean.  Three  types  of  particles  have  been 
observed.  Type  one  particles  are  found  in  agglomerates  of  nominally  30  small 
particles,  type  two  in  agglomerates  of  nominally  5  intermediate  size 
particles  and  type  three  particles  are  completely  dlspersable  singlets.  The 
diameters  of  type  one  and  two  agglomerates  are  approximately  equal  to  the 
largest  type  three  particles. 

V.4.C.  Zeta  Potential  Measurements 

The  apparent  ^-potential  must  be  interpreted  with  care  because 
agglomeration  affects  the  agglomerate's  total  charge  and  drag  coefficient 
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through  Its  cross-sectional  area.  As  will  be  shown  In  the  next  section, 
dilute  dispersions  of  both  pure  and  boron  doped  SI  powders  did  not  exhibit 
continued  agglomeration  In  n-propanol  when  the  water  content  was  less  than 
275  ppm  by  weight.  Thus,  the  calculated  5-potentials  provide  a  qualitative 
measure  of  surface  potential  even  if  the  population  contains  agglomerates. 

The  results  in  Table  V.3  show  that  water  causes  5  to  Increase  for  both 
pure  and  boron  doped  powders  and  that  5  for  pure  SI  Is  constantly  higher  than 
for  the  boron  doped  Si.  The  latter  result  is  surprising  because  FTIR  and 
wetting  experiments  showed  that  the  boron  doped  Si  has  a  relatively  high 
affinity  for  water.  The  effect  of  adsorbed  water  at  the  particle-liquid 
interface  on  5-potentlal  can  be  interpreted  as  resulting  from  the  partial 
hydroxylation  of  the  surface  followed  by  partial  deprotonation.  This  results 
in  a  negative  surface  charge  as  is  the  case  for  silica  in  aqueous 
dispersions.  Hunter17  reports  a  5-potential  of  -110  mV  for  silica  at  low 
ionic  concentration  (10-i*m).  The  Si  5-potential  remains  below  this  value 
because  it  is  expected  that  hydroxylation  will  remain  incomplete  throughout 
the  water  concentration  range  investigated.  The  trend  of  increasing 
5-potential  with  increased  water  concentration  is  consistent  with  these 
expectations.  It  is  concluded  from  these  results  that  the  dominant  charge 
generating  mechanism  in  this  system  is  a  surface  chemical  reaction. 

V.4.D.  Dispersion  Stability 

The  sedimentation  of  dispersions  and  the  agglomeration  rate  of  samples 
prepared  as  described  in  Section  V.3.B.  were  observed.  Weight 
concentrations  of  silicon  (227-SB)  in  n-propanol  ranged  from  4  ppm  to  3X  in 
the  stock  n-propanol  solution.  With  higher  water  contents,  the  Si 
concentration  was  maintained  at  4  ppm. 

The  sedimentation  rates  were  measured  for  stock  n-propanol  containing 
275  ppm  HjO.  A  sediment  layer  was  observed  after  24  hours  for  the  highest 
concentration  of  silicon;  a  clear  supernate  developed  with  longer  times. 

Lower  concentrations  of  silicon  produced  sediment  at  longer  times;  a  fine 
layer  appeared  in  the  10  ppm  dispersion  after  1  week.  Clear  supernate  did 
not  form  in  the  suspensions  when  the  initial  silicon  concentrations  were  less 

than  0.1  wt.Z,  instead  a  cloudy  suspension  remained  for  months.  Quantifying 

35,36 

previous  observations,  high  stability  was  demonstrated  at  low 

concentrations  but  not  at  high  concentrations. 


TABLE  V.3 
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ZETA  POTENTIALS  07  SILICON  POWERS 
IN  N-PROPANOL  WITH  VARIED  WATER  CONCENTRATION 


Silicon  Sample 
(Run  if) 

Wa  ter 
(ppm) 

5-P0  tentlal 
(mV) 

226-S 

275 

-60 

226-S 

40 

-33 

227-SB 

275 

-55 

227-SB 

40 

-29 
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.24.  The  effect  of  water  on  the  rate  of  agglomeration  as  measured  by 
photon  correlation.  Sample  227-SB,  4ppm  In  n-propanol. 


Figure  V.24  shows  the  agglomerate  diameters  in  boron  doped  Si  powder 
227SB  measured  as  a  function  of  time  by  photon  correlation  for  three  water 
contents.  The  rate  of  agglomeration  increases  with  water  content.  In 
agreement  with  the  sedimentation  results,  a  4  ppm  Si  sample  in  the  stock  n- 
propanol  exhibited  no  agglomeration  over  the  period  examined;  samples  of  this 
type  did  not  sediment  for  periods  in  excess  of  1  week.  Pure  Si  powder,  645S, 
also  did  not  exhibit  agglomeration  in  the  stock  n-propanol.  The  enhancement 
of  agglomeration  rate  with  increased  water  content  has  been  reported. 37-39 

The  above  phenomena  can  be  explained  by  considering  certain  implications 
of  the  DLVO  model.  First  it  was  shown  in  Section  V.2  that  the  double-layer 
thickness  for  pure  n-propanol  is  large,  nominally  5  ym.  The  average 
Interparticle  distance  is  800A  for  500A  diameter  Si  particles  at  a  3  wt.% 
concentration.  This  interparticle  spacing  is  well  within  the  double-layer 
thickness  and  is  near  the  maximum  potentials  indicated  in  Figures  V.3  through 
V.6.  These  maximum  potentials  provide  no  effective  barrier  to  agglomeration, 
so  the  agglomerates  probably  grow  to  large  proportions,  settle  out  of 
suspension,  thus  leaving  a  clear  supemate. 

For  lower  solids  concentrations,  a  small  potential  barrier  exists  but  is 
easily  overcome  by  thermal  agitation.  Many  agglomerates  grow  to  a  large  size 
and  settle  out.  Some  however  remain  small  and  remain  suspended  for  extended 
periods  resulting  In  turbid  supernates.  Even  with  very  dilute  dispersions 
and  large  interparticle  separations,  the  diffuse  potential  barrier  Is  not 
great  enough  to  restrict  particle  approach.  The  stability  results  from 
another  factor. 

The  collision  frequency  between  gas  molecules  (Z)  is  given  by: 

Z  -  1/2  n  itT2  c  n2 

where  t  ■  collision  diameter;  c  -  mean  velocity  and  n  -  molecular 
concentration.  The  same  relation  holds  for  dispersed  particles  moving  under 
the  Influence  of  Brownian  motion.  The  rate  of  agglomeration  is  proportional 
to  the  collision  frequency  between  particles;  or  is  proportional  to  n2. 

Dilute  dispersions  can  be  stable  while  concentrated  dispersions  are  not 
simply  because  the  collision  frequency  is  low.  If  small  agglomerates  remain 
after  substantial  sedimentation  has  occured,  the  dilute  supemated  dispersion 
may  appear  stable.  However,  in  initially  concentrated  dispersions,  the 
agglomerates  grow  rapidly  to  very  large  size.  Because  no  substantial 
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fraction  of  small  agglomerates  will  exist  a  clear  supemate  remains  upon 
sedimentation. 

V.5.  CONCLUSIONS 

The  DLVO  model  examined  In  this  paper  predicts  the  unstable  behavior  of 
silicon  In  n-propanol  at  high  concentrations.  The  double  layer  remains  very 
diffuse  due  to  the  low  concentration  of  counter-Ions  in  the  medium  and  the 
potential  barrier  necessary  for  stability  is  not  established.  The 
manipulation  of  surface  charge  or  particle  size  within  reasonable  limits  does 
not  result  in  stability.  For  this  reason  sterlc  means  of  stabilization  Is 
recommended. 

The  application  of  the  Void  method  of  calculating  van  der  Waals 
attractive  potential  for  silicon  particles  with  adsorbed  water  In  n-propanol 
predicted  reduced  stability  as  compared  with  the  non-water  bearing  system. 
These  results  are  in  accord  with  experimental  findings  using  photon 
correlation  and  TEM. 

Photon  correlation  (PC)  was  compared  with  TEM  as  a  means  of  determining 
agglomerate  size  distribution.  Good  agreement  was  found  establishing  PC  as  a 
rapid  method  of  determining  this  parameter.  Application  of  this  method  was 
demonstrated  in  determining  the  effect  of  water  on  dispersion  stability. 
Agglomerate  size  distribution  was  found  to  be  log-normal  by  TEM. 

FTIR  spectra  of  boron  and  non-boron  doped  silicon  powders  were  examined. 
Both  samples  contained  lattice  oxygen  and  SiH.  Physically  adsorbed  water  was 
found  on  the  surface  of  boron  doped  powder,  but  not  on  the  pure  Si  powder. 
FTIR  was  shown  to  be  a  rapid  means  of  chemical  characterization  of  silicon 
powders. 

Microelectrophoresis  was  applied  to  dispersions  of  both  doped  and 
undoped  powders  with  both  moderate  and  low  concentrations  of  water.  Doping 
had  little  effect  on  mobility,  whereas  mobility  Increased  sharply  with 
increased  water. 
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CHAPTER  VI 


SURFACE  TENSIONS  OF  ALUMINA-CONTAINING  LIQUIDS 
J.  M.  Llhraann  and  J.  S.  Haggerty 

ABSTRACT 

The  surface  tensions  of  CC>2  laser-melted  alumina-containing  liquids  have 
been  measured  by  a  static,  pendant-drop  technique  In  air,  helium,  and  a  90% 
helium  and  10Z  hydrogen  mixture.  The  influences  of  up  to  10  wt.%  Cr203,  MgO 
and  TiO^  have  been  studied.  Surface  tension  data  are  also  reported 
throughout  the  whole  concentration  range  between  liquid  alumina  and  zlrconia 
in  air.  An  improved  analytical  technique  based  on  comparing  calculated  and 
observed  melt  contours  Is  reported. 

VI. 1.  INTRODUCTION 

Except  for  salts  and  glasses,  relatively  few  data1”5  exist  for  the 
physical  and  chemical  properties  of  liquid  ceramics.  This  paper  reports 
static  surface  tensions  of  alumina-containing  melts  measured  by  the  pendant 
drop  technique.  Future  dynamic  techniques  will  provide  independent 
measurements  of  the  surface  tension-to-densi ty  ratios  as  well  as  the  melt 
viscosities.  The  Al203~Cr203,  AljOj-MgO,  AIjOj-TIOj  and  AljOj-ZrOj  systems 
were  selected  because  at  present  they  represent  the  commercially  most 
important  and  technically  most  interesting  alumina-based  ceramics. 

VI. 2.  ANALYTICAL 

Broken  and  distorted  bonds  between  atoms  at  the  llquld-vapor  Interface 
cause  an  excess  energy  defined  as  the  surface  energy.  For  incompressible 
liquids,  the  surface  energy  Is  equal  to  the  surface  tension  because  surface 
stresses  are  relieved  by  viscous  flow.  The  surface  tension  causes  a  pressure 
difference,  AP,  across  a  curved  Interface,6”7  defined  by: 

AP  -  yd/RjU)  +  1/R2<*>> 
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(1) 
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where  y  Is  the  surface  energy  or  tension  of  the  liquid-vapor  interface,  and 
Rj  and  R2  are  the  principle  radii  of  curvature  of  the  interface  at  an 
altitude  z.  The  internal  pressure  within  a  fluid  body  having  a  density  p 
varies  with  altitude  as  pgz;  the  changing  pressure  causes  the  local  shape  to 
vary  accordingly. 

The  drop  shape  is  used  to  define  the  surface  tension  in  the  pendant-drop 
technique.  Figure  VI. 1  shows  the  pendant-drop  geometry  and  Equation  2 
defines  the  analytical  relationship  between  the  liquid  properties  and 
pendant-drop  shape: 

2y/b  -  y(l/Rj (z)  +  1/R2<z))  =  pgz  (2) 

where  b  Is  the  drop  radius  at  z«0. 

Pendant-drop  shapes  can  be  analyzed  using  the  differential  form  of 
Equation  2: 


(£i^i)(£) 
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f - thllll - )  +  f - 2z/dr - ^  (3)  e 

(1  +  (dz/dr)2)3/2  r(l  +  (dz/dr)2 )* /2  . 


Solutions  to  Equation  3  are  families  of  curves  parameterized  by  6 
(Figure  VI. 1)  where: 

6  *  -pgb2/y  (4) 

and  r/b  and  z/b  values  are  tabulated®-^  as  a  function  of  8.  Although 
Bashforth  and  Adams'  tables8  contain  exact  solutions  to  Equation  3,  they  have 
not  been  extensively  used  since  they  require  the  experimental  determination 


of  b. 


Alternatively,  the  pendant  drops  can  be  analysed1 1-1 3  using  the 


relationship 


where  H  is  a  numerical  factor  tabulated  as  a  function  of  S  *  r  /r 

s  e 

(Figure  VI. 1).  This  procedure  is  simpler  and  potentially  more  accurate  since 
It  eliminates  the  necessity  of  measuring  b.  The  method  is  still  subject  to 
considerable  error  since  it  depends  on  the  rg  value  measured  at  an  altitude 


Stauffer14  analyzed  the  errors  expected  using  this  analytical  procedure. 
Assuming  a  1%  measurement  error,  he  showed  that  the  uncertainty  in  y/p  is  20% 
for  S  *  0.4  while  It  decreases  to  2.6%  for  S  ■  0.85.  The  extreme  sensitivity 
of  the  technique  to  measurement  errors  is  evident  from  the  small  differences 
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In  the  positions  of  curves  A,  B  and  C  in  Figure  VI. 1,  at  an  altitude 

z  m  2r  • 
e 

The  use  of  longer  length-to-d lame  ter  ratios  improves  the  precision 

because  the  pendant-drop  contours  separate  from  one  another  as  shown  in 

Figure  VI. 1.  The  accuracy  also  improves  because  the  actual  pendant-drop 

shape  is  less  subject  to  distortions  caused  by  non-circular  cross-sections, 

non-planar  solid-liquid  Interfaces  and  non-freely  hinging  solid-liquid 

interfaces.  Many  authors*2-15  have  recommended  using  pendant-drop  lengths 

that  exceed  3. 0-4.0  r  . 

e 

We  found  this  procedure  gave  unacceptable  accuracy  with  the  oxides 
investigated.  In  fact,  pendant-drop  lengths  in  excess  of  2. 2-2. 6  r^  could 
not  be  achieved  with  the  y/p  values  characteristic  of  this  work,  because 
Rayleigh  instability  caused  the  pendant-drops  to  neck  off  when  drop  lengtb- 
to-rsratios  exceeded  a  critical  value  («=  3.0)  apparantly  related  to  the  zone- 
height  Instabilities16  found  for  floating-zone  crystal  growth  between  unequal 
rod  diameters.  Longer  drop  lengths  would  require  lower  y/p  values.  Also,  we 
were  not  able  to  reliably  measure  rg  values  with  1.02  precision  because  of 
the  combined  effects  of  r  errors,  absence  of  a  sharp  image  boundary  caused 
by  optical  flare  and  an  uncertain  position  caused  by  vibrations  that  appeared 
Intrinsic  to  the  high- temperature  melts. 

Therefore,  we  elected  to  compare  shapes  computed  from  Equation  3 
directly  with  observed  shapes  to  minimize  the  errors  resulting  from  short 
pendant-drops.  By  making  a  visual  comparison  over  the  entire  lengths  of  the 
contours,  it  was  possible  to  make  a  precise  selection  of  the  correct  computed 
curve. 

The  procedures  used  to  calculate  the  pendant-drop  shapes  followed  those 
given  in  References  17-19.  The  entire  contour  Is  defined  by  the  assumed  rg 
and  6  values;  8  is  defined  by  rg/b.  For  an  assumed  rfi  value,  a  family  of 
contours  was  generated  as  a  function  of  b.  Typical  such  contours,  A-H,  are 
shown  in  Figure  VI. 1  for  a  particular  value  of  rfi.  Iteration  over  a  range  of 
revalues  provided  a  complete  set  of  curves,  each  of  which  correspond  to  a 
specific  b2/8.  Equation  4  is  used  to  calculate  y/p  for  the  best  fitting 
curve. 

It  is  important  that  the  observed  and  calculated  contours  be  the  same 

over  their  entire  lengths.  Figure  VI. 1  illustrates  the  expected  shapes  for 

constant  r  values  with  different  y/p  values.  Curves  A-C,  Figure  VI. 2, 
e 
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Figure  VI. 1.  Shows  some  pendant-drop  shapes  computed  for  a  constant  r£ 

value.  Curves  D  through  H  in  this  order  correspond  to  (3  values 
(Equation  3)  of  -0.25,  -0.31,  -0.37,  -0.475,  and  -0.55.  Curves 
A,  B  and  C,  though  defined  by  very  close  values  of  S  -r  /rp 
(Equation  6)  and  almost  undlstlnguishable  at  the  altitude 
z  ■  2re,  correspond  to  y/p  values  as  different  as  214.73, 

201.79  and  187.24,  respectively. 


Figure  VI. 2.  Illustrates  the  non-linear  aspect  of  Equation  3.  Parallel 
curves  like  A,  B,  and  C  correspond  to  very  different  y/p 
values:  190.4  (A),  253,7  (b)  and  352.5  (C),  whereas  non- 
parallel  contours  D,  E  and  F  define  similar  values  of  v/o: 
249.5  (D)  253.7  (E),  258.2  (F). 
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Illustrate  the  error  In  y/p  resulting  from  selecting  a  contour  that  parallels 
the  observed  pendant-drop  shape.  Curves  D-F,  Figure  2,  Illustrate  that  the 
pendant-drop  shapes  are  In  fact  quite  a  strong  function  of  rg  for  similar  y/p 
values.  With  a  little  practice,  it  Is  possible  to  distinguish  between 

matching  and  incorrect  contours. 

20  21  5 

Except  for  pure  AljOj  ’  and  Cr203  ,  the  densities  of  the  investigated 
liquids  have  not  been  reported.  Rule  of  mixture  approximations  and 
assumptions  given  in  Reference  3  were  used  to  calculate  liquid  densities  of 
Al203-Cr203,  Al203-Mg0,  Al203-T102  and  Al203~Zr02  mixtures.  Our  calculations 
assumed  the  density  of  liquid  AljOj  at  the  melting  point  to  be  3.01  g/cm3. 

The  densities  of  other  end-member  molten  oxides  assumed  a  20. 4%  volumetric 
expansion  upon  melting.  Summarized  in  Figure  VI. 3  are  the  p  results  used  in 
y  calculations. 

VI. 3.  EXPERIMENTAL 

Experimentally,  melts  were  created  on  the  bottom  ends  of  polycrystalline 
feed  rods  in  a  controlled  atmosphere  chamber.  The  incandescent  melts  were 


Cone. of  additive  in  wt.% 


Figure  VI. 3.  Density  of  the  various  melts  as  calculated  with  the 
approximations  detailed  in  Section  VI. 2. 
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projected  with  a  single  lense  (20x  magnification)  onto  a  surface  where  they 
were  compared  with  calculated  contours. 

Heating  was  accomplished  with  four  10.6  pm  wavelength  beams  of  light 
that  orthogonally  impinged  onto  the  melts  in  a  radial  direction.  The  light 
source  was  a  1350  watt-C02  laser  custom-designed  to  emit  two  beams.  Each  of 
the  emitted  beams  was  divided  into  two  beams  with  roof-prisms.  The  beam 
dimensions  on  the  melt  surfaces  were  controlled  by  the  relative  positions  of 
the  focusing  lenses  and  the  melt.  Typical  spot  diameters  equaled  the  feed- 
rod  diameters. 

Laser  heating  is  particularly  applicable  to  high  temperature  materials. 
The  laser  has  no  characteristic  temperature  to  set  maximum  temperature 
limits.  Also  ambient  atmospheres  can  be  selected  without  serious 
restrictions. 

Feed  rods  were  isos ta tically  pressed  from  powders  batched  to  the  desired 
overall  compositions  with  mixtures  of  pure  oxide  powders.  Nominally  3  mr 
diameter  rods  were  sintered  for  12  hours  at  1200°C,  typically  producing  a  65'= 
dense,  10  pm  grain  size  rod.  Further  densification  occurred  over  a  distance 
of  nominally  one  rod  diameter  when  they  were  melted.  No  material  was 
withdrawn  from  the  melt  onto  the  feed  rod,  so  the  melt  composition  was  equal 
to  the  batch  composition.  With  the  exception  of  the  high  ZrC>2  compositions, 
vaporization  losses  were  negligible. 

The  relative  beam  positions  on  the  melt  and  beam  powers  were  adjusted  to 
make  the  solid-liquid  interface  as  flat  as  possible.  When  properly  adjusted, 
the  altitude  of  the  solid-liquid  interface  was  uniform  within  ±  0.05  mm. 

Shaft  rotation  flattened  the  interface  further  but  could  not  be  used 
effectively  in  surface  tension  measurements  because  vibrations  made  the  drop 
shape  uncertain  and  severely  limited  the  maximum  drop  length  that  could  be 

achieved  before  it  fell  from  the  feed  rod. 

22  23 

Some  superheating  ’  occurred  where  the  laser  beams  intersected  the 
melts.  By  expanding  the  beam  diameters  to  equal  the  feed  diameters,  the 
radiant  flux  onto  the  melt  surfaces  was  made  as  uniform  as  possible.  This 
adjustment  combined  with  the  action  of  a  spherical  radiation  shield 
positioned  to  refocus  reflected  and  emitted  radiation  on  the  center  of  the 
melt  minimized  the  superheat  to  *  20°C.  Free  convection  in  the  melt  acted  to 
reduce  temperature  gradients  and  insure  compositional  uniformity. 


i  -  . 
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Melt  contours  were  analysed  with  progressively  Increasing  melt  volumes 
until  the  melt  dropped  off  to  maximize  the  length-to-diameter  ratio.  The 
final  S  values  were  *  0.85.  For  this  study,  300  theoretical  contours  were 
computed  for  15  values  of  rg  ranging  from  2.17  to  A. 07  mm  and  20  values  of 
(-6)  ranging  from  0.25  to  0.55.  The  incremental  y/p  between  successive 
curves  is  approximately  1%.  Generally  each  data  point  represents  an  average 
of  5  measurements. 

VI. 4.  RESULTS  AND  DISCUSSION 

The  results  demonstrate  the  superior  precision  of  this  technique  for 

analyzing  pendant-drops.  The  extreme  between  y/p  values  for  5  independent 

measurements  was  typically  1 7,  when  determined  by  comparing  contours.  This 

reproducibility  is  comparable  to  the  incremental  differences  between  the 

calculated  contours  used  to  analyse  the  experimental  contours.  In  contrast, 

y/p  values  determined  from  measured  r  and  r  values  varied  as  much  as  ±10%. 

e  s 

The  accuracy  potential  provided  by  this  technique  exceeds  our  knowledge 

of  the  melt  densities.  For  pure  alumina,  reported  liquid  densities  at  the 

24  25 

melting  point  range  from  2.51  to  3.05  g/cm3  .  We  adopted  the  more  recent 

_  1  27 

and  accepted  value  of  3.01  g/cm3.  ’  Since  the  densities  of  other  studied 
liquids  have  not  been  characterized  experimentally,  we  calculated  them  using 
an  average  value  of  the  thermal  expansion  coefficient  of  the  solid  between 
the  ambient  temperature  and  the  melting  point,  and  assuming  a  20.4%  volume 
expansion  upon  melting.  It  is  evident  that  the  y-accuracy  is  limited  by  the 
p  accuracy  at  this  time.  We  would  speculate  that  our  y-errors  are  in  the 
range  of  2-3%  (13-20  ergs/cm2).  As  can  be  seen  from  Table  VI. 1,  this 
precision  is  much  better  than  those  obtained  so  far  with  other  surface 
tension-measurement  technique.  The  difference  between  our  measured  value  for 
the  surface  tension  of  A1203  In  He  and  that  cited  in  Reference  5  (570 
ergs/cm2)  exceeds  this  estimate  of  error.  This  difference  may  reflect  a 
systematic  difference  between  the  pendant  drop  and  bubble  pressure 
measurement  techniques. 

The  results  of  y  determinations  for  Al203-Cr203,  AljOj-MgO,  Al203~Ti02 
and  Al203~Zr02  systems  for  each  of  the  three  atmospheres  investigated  are 
shown  in  Figures  VI. 4  through  VI. 7.  The  y-axes  are  truncated  to  expand  the 


113 


TABLE  VI. 1 

CITED  SURFACE  TENSION  VALUES  OF  LIQUID  ALUMINA 


Surface  Tension 
(ergs/cm2 ) 


665± 15 
625± 14 
610± 13 


577+80 

586+81 


690+48 


574+68 


Atmosphere 


air 

He 

He+10%H, 


vacuum 

vacuum 


not  specified 


vacuum 


Method  Used 

Reference 

pendant  drop 
(p=3. 01 ) 

this  study 

theoretical  value 
(based  on  the  (electrostatic) 
forces  between  ions  in  a 
coordination  or  nearly  ordered 
mel  t) 

Sokolov2  £ 

drop-weight  method  (p=2.5) 

drop-weight  method  (p=3.2) 

Von 

War  tenberg2** 

pendant  drop  (p=2.97) 

Kingery2 

sessible  drop,  from 
contact  angle  value 

Bartlett2  5 

pendant  drop  (p=3.01) 

Rasmussen2  7 

maximum  bubble  pressure 
(p-3.03) 

Anisimov5 
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Figure  VI. 4.  Experimental  surface  tension  of  the  liquids  in  air. 


Figure  VI. 5.  Experimental  surface  tension  of  the  liquids  in  helium. 
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ATMOSPHERE : 
90%He  - 10%  H2 


j _ L 


to 


Pure  AbO^  Cone,  of  nrlrlitivp  in  wt  % 


Figure  VI. 6.  Experimental  surface  tension  of  the  Liquids  in  a  90°.  heLium 
10%  hydrogen  mixture. 


Figure  VI. 7.  Experimental  surface  tension  of  various  Al203-Zr02  liquids  in 


scales.  The  Al2C>3-Zr02  data  were  plotted  separately  front  the  other  data 
taken  In  an  air-atmosphere  because  ranges  of  variables  were  wider. 

In  general,  surface  energies  increased  with  increasing  oxygen  partial 
pressure.  For  pure  AljO^ ,  y  was  665,  625  and  610  ergs/cm2  for  air.  He  and 
He-H2  atmospheres  respectively. 

The  effect  on  surface  energies  of  additions  up  to  ~  10  weight  percent 
was  generally  either  a  small  monotonic  increase  or  decrease.  Ti02,  the 
exception,  caused  a  rapid  decrease  in  y  for  levels  <0.5  wt.Z,  followed  by 
essentially  constant  y  values  for  higher  concentrations.  All  samples 
characterized  by  electron  microprobe  showed  identical  surface  and  bulk 
compositions. 

For  an  air  atmosphere,  additives  generally  caused  the  y  values  to 
increase  or  decrease  in  the  same  direction  that  they  cause  the  liquidus 
curves  to  change.  For  small  concentrations,  C^203  ls  t*ie  on^y  additive  that 
raises  the  melting  point  and  is  the  only  additive  that  raises  the  y  value. 

The  minimum  y  value  for  the  Al203~2r02  system  corresponds  to  the  eutectic 
composition,  approximately  42  wt.%.  For  Zr02  contents  >  70  wt.%,  y  decreases 
with  increasing  Zr02 ,  violating  this  generality.  This  departure  may  result 
from  a  thermal  reduction  of  the  high  ZrOj-content  melts  as  evidenced  by  the 
darkening  of  the  quenched  pendant-drops,  following  the  trend  of  reduced  •, 's 
with  reduced  oxygen  partial  pressures. 

The  effects  of  additives  on  y's  in  He  and  He-H2  atmospheres  are  more 
complex;  their  magnitudes  do  not  follow  trends  expected  on  the  basis  of 
oxygen  partial  pressure.  No  data  was  found  for  the  effect  of  oxygen  partial 
pressure  on  melting  points,  so  this  correlation  could  not  be  examined.  It  is 
apparent  that  y's  for  these  systems  are  sensitive  to  atmosphere  and 
composition  even  though  the  causes  are  not  understood. 

VI. 5.  SUMMARY  AND  CONCLUSIONS 

An  improved  technique  for  analyzing  the  surface  energies  of  liquids  in  a 
pendant-drop  configuration  has  been  developed.  The  procedure,  based  on 
comparing  observed  and  calculated  contours,  is  particularly  applicable  to 
melts  that  cannot  be  extended  to  high  length-to-d lame  ter  ratios  and  to  high 
temperature  melts  for  which  precise  dimensional  characterization  ls 
difficult.  The  results  demonstrated  the  technique  with  several  alumina 
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containing  melts.  The  precision  of  the  y/p  measurements  appears  to  be 
approximately  1%.  The  accuracy  of  the  y  values  is  limited  by  the  lack  of  p 
values. 

Both  composition  and  atmosphere  have  strong  effects  on  y  values.  In  an 
air  atmosphere,  the  effect  of  additives  on  y'e  parallels  their  effect  on 
melting  points.  In  He  and  He-H2  atmospheres  the  effects  are  more  complex. 
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